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ABSTRACT 
The  direct  problem  in  transient  heat  conduction  requires  the 
determination  of  conditions  at  an  interior  location  when  conditions 
are  known  at  the  boundaries  of  a  solid.   Conversely,  the  inverse 
problem  requires  the  determination  of  conditions  at  the  boundaries 
of  a  solid  when  conditions  are  known  at  an  interior  location.   Con- 
sequently special  methods  are  required  in  the  solution  of  the  inverse 
problem.   A  new  method,  numerical  inversion  of  the  Laplace  transform, 
is  used  to  solve  this  complex  problem.   Application  of  this  numerical 
technique  to  the  semi-inf init^pWHBP*>Wftg"  cylinder,  and  sphere  is 
made,  and  the  accuracy  of  solution  is  discussed.   This  method  of 
solution  provides  the  engineer  with  a  simple,  powerful  tool  that  can 
be  used  in  the  determination  of  heat  transfer  phenomena  in  a  solid. 


Thesis  by:   Terrill  Jay  Wendt,  entitled  "Application  of  Numerical  Inversion 
of  the  LaPlace  Transform  to  the  Inverse  Problem  in  Transient  Heat  Conduction". 
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I .     INTRODUCTION 

Traditionally,  heat-transfer  problems  are  concerned  with  the 
determination  of  conditions  at  an  interior  location  when  conditions 
are  known  at  the  boundaries  of  a  solid.   Such  problems  are  known  as 
direct  problems,  and  their  solutions  are  obtained  through  straight- 
forward  classical  techniques  [5].   In  many  applications  (i.e., 
re-entry  thermodynamics,  rocket  motors,  quenching),  it  is  desired 
to  determine  the  conditions  at  the  boundaries  of  a  solid  when 
conditions  are  known  at  an  interior  location.   Problems  of  this  type 
are  known  as  inverse  problems  and  require  special  methods  of  solution. 

Many  solutions  to  the  inverse  problem  of  heat  conduction  exist 
in  the  literature.   Stolz  [19]  uses  numerical  inversion  of  a  specific 
direct  problem  to  solve  for  the  heat  flux  at  the  surface  of  a  solid. 
He  notes  that  there  is  appreciable  damping  or  lag  in  time  of  the  heat 
flux  from  the  interior  location  to  the  surface.   This  damping  resulted 
in  oscillations  of  the  approximate  solution  with  resultant  errors. 
Burggraf  [4]  obtains  an  "exact"  solution  to  the  inverse  problem  in 
heat  conduction  with  internal  heat  generation.   However  Beck  [4]  in 
commenting  on  this  solution  notes  that  this  "exact"  method  must 
utilize  inexact  measured  temperature  data  for  which  accurate  time 
derivatives  (necessary  for  this  "exact"  solution)  are  very  difficult 
to  obtain.   In  commenting  on  the  inaccuracies  involved  in  using 
measured  temperature  data,  Frank  [12]  suggests  using  a  least  squares 
fit  to  the  data  available  in  order  to  obtain  an  essentially  smooth 
temperature  function  for  use  in  the  solution  of  the  inverse  problem. 
Sparrow,  jst.  al_.  [18]  use  a  modification  of  Laplace  transform  inversion 


* 
Numbers  in  brackets  refer  to  references  in  the  bibliography. 
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by  applying  a  convergence  factor  to  the  temperature  function  in  the 
s  plane.   Deverall  and  Channapragada  [11]  utilize  a  convergence  factor 
to  obtain  a  solution  for  the  heat  flux  in  an  inverse  problem. 

It  is  intended  here  to  present  a  simple  straightforward  method 
of  obtaining  solutions  to  inverse  problems  of  practical  interest.   The 
numerical  technique  involved  in  inverting  the  Laplace  transform  (used 
in  solving  the  inverse  problems)  was  developed  by  Bellman  and  Kalaba 
[2  and  3].   This  method,  based  on  orthogonal  function  theory  and 
Gaussian  quadrature,  is  developed  and  tested  for  accuracy  (Tables  I 
and  II) .   This  numerical  technique  is  applied  to  the  semi- infinite 
solid,  the  "long"  cylinder  and  the  sphere.   The  accuracy  of  solution 
is  compared  graphically  with  the  known  exact  solutions.   The  FORTRAN  IV 
programs  used  in  the  solutions  are  listed  in  Appendix  A  for  use  by 
the  reader. 
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II.    NUMERICAL  TECHNIQUE 

INTRODUCTION 

The  numerical  technique  Involved  was  developed  by  Bellman  and 
Kalaba  [2]  and  involves  the  application  of  orthogonal  function  theory 
to  the  inversion  of  the  Laplace  Transform  as  opposed  to  using  the 
Convolution  Theorem  which  involves  integration  in  the  complex  plane. 

It  is  important  to  note  that  the  method  developed  in  this 
section  can  be  applied  to  other  problems  in  physics,  biology,  eco- 
nomics, chemotherapy  and  radiation  [2]. 

THE  LAPLACE  TRANSFORM 

We  will  first  consider  some  of  the  properties  of  the  Laplace 
transform  and  limitations  involved  in  applying  the  inverse. 

Let  f(t)  be  any  sectionally  continuous  function  defined  for 
t  ^>  0.   Then  the  Laplace  transform  of  f (t)  is  defined  as 

F(s)  =  L{f(t)}  = 


e  St  f(t)dt  (1) 


o 
where  the  real  part  of  s  is  large  enough  to  make  the  integral  (1) 

convergent. 

The  Laplacian  Operator  L  is  a  linear  operator  and  possesses 
many  desirable  properties,  mathematically  and  computationally.   Some 
of  these  properties  are  illustrated  in  the  following  theorems  which 
are  included  here  for  completeness: 

THEOREM  I:    If  L{f  (t)}  =  F  (s)  and  L{f  (t)}  =  F  (s) 

11  2  2 

then  L{f  (t)  +  f  (t)}  =  F  (s)  +  F  (s) . 

12  12 

THEOREM  II:    If  L{f(t)}  -  F(s) 

then  Lff^-}  =  sF(s)  -  f(0). 
d  t 
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THEOREM  III:    If  L{f(t)}  =  F(s) 


then  LI 


9Nf(t)[   d1 


N 


dx 


N 


F(s) 


THEOREM  IV:    If  L{f(t)}  =  F(s) 


c 


f(t)dt >   =  -  F(s) 


THEOREM  V: 


If  L{f(t)}  =  F(s) 


then  L{f(at)}  =  -  f(—). 
a   va/ 

THEOREM  VI:    Initial  Value  Theorem. 

If  L{f(t)}  =  F(s) 


then 
Limit 


_,,  ,    Limit  -,.% 
sF(s)  =       f(t) 


THEOREM  VII:    Final  Value  Theorem. 
If  L{f(t)}  =  F(s) 
then 

Limit  sF(s)  =  Limit  f(t). 

g->oo       v  '       fc-XJO 

THEOREM  VIII:   Superposition  Theorem. 

If  L{f(t)}  =  F(s)    and   L{g(t)}  =  G(s) 
then 


rt 


f(t)g(t-t)dt?  =  L 


g(t)f(t-t)dt 


«  FCs)G(s) 
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INSTABILITY  OF  THE  INVERSE  OF  THE  LAPLACE  TRANSFORM 

The  instability  of  the  inverse  of  the  Laplace  transform  can 
be  illustrated  by  considering  the  function 

f(t)  =  sin  cot. 
The  Laplace  transform  of  f (t)  is  then 

L{sin  ait}  =  -^2  (2) 

Now  it  is  known  that  sin  0Jt  is  bounded  in  amplitude  by  (1,-1);  but  as 
co  is  increased,  sin  cot  oscillates  correspondingly  more  rapidly. 
However  the  Laplace  transform  of  sin  cot  approaches  zero  as  co  is 
increased. 

Limit   co 


s2+co2 

co-*» 


=  0 


In  other  words  small  changes  in  s  can  cause  arbitrarily  large 
changes  in  t  with  respect  to  the  period  of  oscillation.   Consequently 
we  cannot  sift  or  filter  out  rapid  oscillations  by  using  values  of 
F(s).   Such  difficulties  as  mentioned  above  can  be  avoided  only  if 
f(t)  is  an  essentially  smooth  function,  which  is  an  implicit  assump- 
tion in  the  following  development. 

NUMERICAL  QUADRATURE 
There  are  several  methods  that  can  be  used  to  approximate  a 
definite  integral  in  terms  of  a  finite  sum.   The  method  employed  here 
is  based  on  Gaussian  Quadrature  and  is  used  because  of  convergence 
properties  and  exactness  of  solution  [15].   Consider  the  integral 

approximation 

rb  M 


w(x)f(x)dx  =  2_  A.f(x.).  (3) 

a  i-1  X   1 


By  minimizing  the  function 


w(x)f(x)  -  2_   A  f(x  ) 
i=l 


(A) 
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we  can  assume  that  the  summation  involved  in  equation  (3)  is  exact 
for  a  prescribed  accuracy,  e,  for  any  polynomial  of  degree  less  than 
or  equal  to  M.   The  x.'s  and  the  A  's  are  to  be  determined  so  that 
equation  (3)  is  exact  for  a  polynomial  f(x)  of  high  degree.  Here 
the  x.'s  are  called  the  roots  (or  nodes),  the  A  's  are  called  the 
coefficients  (or  weights) ,  and  w(x)  is  called  the  weighting  function 
of  the  function  f(x). 

With  this  background  let  us  examine  the  numerical  approxima- 
tion to  an  integral  which  will  be  used  later  to  invert  the  Laplace 

transform.   Let 

M 

(5) 


7T(x)  =  1   "  (X-X  )  , 

i=l 
and  let  PM(x)  be  a  polynomial  of  degree  N.   Then  we  have 


PN(x)  =  TT(x)QN_M(x)  +  PM-1(x). 
Since   PN(x.)  =  P^U.) 


(6) 


i  =  1,2 M 


then 


w(x)PN(x)dx  = 


fb 


w(x)PM_1(x)dx  + 


w(x)7T(x)QN_M(x)dx.  (7) 


Now  let  us  choose  the  x.'s  such  that 

l 


w(x)TT(x)QN_M(x)dx  =  0, 


(8) 


for  all  polynomials  QN_M(X) •   Then  from  equation  (3)  we  have 


w(x)PN(x)dx  =  2_   A.PN(Xi) 


(9) 


a  i=l 

If  equation  (8)  is  to  be  valid  for  all  polynomials  Q-T  M  for 


<N-M 


N  _>  M,  it  must  be  valid  for 

rb 

K 
w(x)tt(x)x  dx  =  0, 

a 


K  -  1,2,...,N-M 


Let  us  now  expand  tt(x)  into  a  polynomial 


7T(X)  =2_  C 

i=o 


(10) 


(11) 
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so   that   equation   (10)   becomes 


b  /M 

w(x)[  >      C±x  /x  dx 

a         Vi=o 


i  =  0,1,2,... ,N-M 


(12) 


which  is  a  system  of  M  equations  in  the  c.'s  if  N  =  2M-1. 

The  system  of  equations  in  (12)  has  a  unique  solution  ^f  and 
only  if  the  matrix  of  coefficients  is  nonsingular. 

The  determinant  of  this  matrix  is  called  a  Vandermonde  deter- 
minant and  can  be  easily  evaluated,  since  w(x)  >  0  in  the  inverval 
(a,b),  as 

1   x   . 


N 


1   *N  ' 


.  x. 


N 


N 


*H 


i>j 


(x.-Xj). 


(13) 


Since  the  x  ' s  are  real  and  distinct  points  (roots),  the  determinant 
in  (13)  does  not  vanish  and  the  b  's  can  be  uniquely  determined. 

It  can  be  shown  that  the  zeros  of  tt(x)  ,  determined  by  its 
coefficients  (c.),  are  real  and  in  the  interval  (a,b).   It  can  also 
be  shown  that  Q   (x)  is  a  polynomial  of  degrees  M-l,  where  N  =  2M-1, 

Now  let  tt.(x)  and  tt„(x)  designate  the  functions  \.(x)  deter- 
mined by  their  coefficients  for  M  -   i  and  M  =  j  respectively.   Then 

rb 


w(x)tt  (x)TT.(x)dx  =  0 


i*  J 


(14) 


If  i  <  j,  then  tt.(x)  is  a  polynomial  of  degree  less  than  j;  and  the 

result  follows  from  equation  (8),  which  has  been  verified.  From 

equation  (14),  77.  (x)  and  tt.(x)  are  said  to  be  orthogonal  on  the 
interval  (a,b)  with  respect  to  the  weight  function  w(x) . 


17 


For  each  value  of  M,  the  x.'s  in  equation  (3)  are  taken  as 

the  zeros  of  tl.(x)  as  determined  by  its  coefficients.   The  A  's  in 
M  i 

equation  (3)  can  then  be  determined  by  the  method  of  undetermined 
coefficients.   It  can  be  shown  that  all  of  the  A.'s  are  positive. 

Therefore  for  a  given  weight  function  w(x)  and  interval  of 
integration  (a,b)  ,  the  zeros  of  the  polynomials  tt  (x)  ,  as  determined 
by  their  coefficients  (c„),  need  be  computed  only  once.   Likewise 
the  corresponding  A.'s  need  be  computed  only  once. 

GAUSS -LEGENDRE  QUADRATURE 
Since  we  have  been  dealing  with  any  orthogonal  function  t  (x)  , 
let  us  indicate  the  specific  orthogonal  functions  used  in  the  fol- 
lowing development.   The  Legendre  polynomials  are  the  unique  Gauss 
quadrature  polynomials  [15]  in  the  interval  (-1,1)  (or  any  finite 
interval)  with  weighting  function 

w(r)  =   1  (15) 

Let  us  denote  the  Legendre  polynomials,  in  accordance  with 

common  practice,  by  the  notation  PN(r)  where 

t,  /  n    /  n   (2N)!  ,,« 

PN(r)  =  7T(r)  2j^(N;)2  ,  (16) 

and  the  interval  of  expansion  is  (-1,1).   These  polynomials  are  defined 

by  the  Rodrigue's  formula 

Vr)  "  W!  J&  (r2-jf  •  (") 

From  equation  (14)  it  follows  that 

'l 


and     PN+1 


PN(r)PM(r)dr  =  0  M  j>  N  .  (18) 

'-l 

The  recursion  relations  for  the  Legendre  polynomials  are  then 

Wr)  ■  (tsi)1  Vr>  -  (m)  h-iM<  (19) 

w-V^fe)'.'"1'  (20) 
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Consequently  the  least  squares  polynomial  approximation 

(Gaussian  Quadrature)  to  f(x)  over  (-1,1)  with  respect  to  constant 

weighting  function  is 
N 
f(r)  =  2>KPK(r),  (21) 

K=0 

1 


u  u  2K+1 

where   b  =  — r — 
K    2 


f(r)P  (r)dr.  (22) 

IS. 

■i 


Therefore  the  coefficients  (or  weights)  A.  are 

dr  ,  (23 


A.  = 

l 


(r-r.)  PMx(r) 

— 1      1    JN 


where  r  '  s  are  the  roots  of  the  Legendre  polynomials  P  (r) . 
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SHIFTED  LEGENDRE  POLYNOMIALS 

We  have  now  examined  a  specific  orthogonal  polynomial,  the 
Legendre  polynomial.   However,  it  is  most  convenient  in  the  following 
development  to  consider  the  shifted  Legendre  polynomials.   The  reason 
for  employing  these  polynomials  is  twofold .   First,  the  interval  of 
expansion  (0,1)  of  the  shifted  Legendre  polynomials  corresponds  to 
that  involved  in  obtaining  an  explicit  inverse  matrix.   Second,  the 
coefficients  of  the  shifted  Legendre  polynomials  are  all  integers, 
which  is  an  important  property  computationally. 

Since  the  roots  and  weights  have  been  tabulated  for  the 
interval  (-1,1)  [19],  it  is  a  simple  matter  to  obtain  the  roots  and 
weights  of  the  shifted  Legendre  polynomials  in  the  interval  (0,1). 

•  Let    * 

PN   (x)  =  PNd-2r)  (24) 

* 
where  P   (x)  are  the  shifted  Legendre  polynomials „   Consequently  the 

roots  of  the  shifted  Legendre  polynomials  are 

1+rj 
x±  -  -yi  ,  (25) 

and  the  weights  are 

A±  -  -|  .  (26) 
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roo 

e"Stu(t)dt.  (27) 


NUMERICAL  INVERSION  OF  THE  LAPLACE  TRANSFORM 

We  have  now  laid  a  proper  foundation  for  obtaining  an  explicit 

inverse  matrix  to  be  used  in  inverting  the  Laplace  transform.  Let  us 

consider  the  approximate  inversion  of  the  Laplace  transform  through 

numerical  techniques. 

Let 

F(s)  = 

o 
Now  we  make  a  change  of  the  variable  of  integration  to 

x  =  e_t  (28) 

therefore  t  =  -  log  x  (29) 

Substituting  equation  (29)  into  equation  (27)  we  obtain 

xS"1u(-log  x)dx  =  F(s).  (30) 

o 
Now  let  a(x)  =  u(-log  x)  (31) 

and   then  equation   (30)   becomes 

xS_:La(x)dx  =  F(s).  (32) 

J 
o 

s-1 
If  we  make  g(x)  =  a(x)x.    ,  then  using  the  quadrature  formula 

'i        J-   *  X 
g(x)dx  =  X  A  g(x  ),  (33) 

i=l 
o 

we  obtain 

N 
5"  A.  x.S  a(x.)  =  F(s).  (34) 

i=l   1   1      X 

Letting  s  =  1,2,...,N,  we  obtain  a  linear  system  of  N  equations  in 

the  N  unknowns,  a(x.),  i  =  1,2,..., N, 

J-    *\ 

>  A.  x  a(x.)  =  F(K+1)     k   =  0,1,..., N-l        (35) 

i=l  X      1 

It  is  our  intention  here  to  obtain 

N-l 
a(x.)  »^V,„F(ic+l),       i  =  1,2 N,        (36) 

1     ^ IK 

K=0 

where  the  a.  's  are  constants  independent  of  F(k+1) ,  determined  once 

IK 


and  for  all  for  each  N,  and  tabulated. 
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Let  us  now  make  another  change  of  variables  such  that 

y±  =  A±*a(Xi),  (37) 

and  a  =  F(k+1) .  (38) 

Equation  (35)  then  becomes 

N 

21  xiKyi  =  \  '  K  -  0,1,.-.,N-1.       (39) 

i=l 

Multiplying  the  K   equation  by  a  parameter  q  ,  and  adding  the 

Is. 

corresponding  terms  in  all  the  equations,  we  get 
N    /N-l     \    N-l 

ZyJ^-Vi  7  =  ^a<V  (40) 

1=1   IK=0      /     K=0 

Now  setting 
N-l 

f(x)  =  ^T<iKxK   ,  (41) 

K=0 

N  N-l 

we  have  ^iy±HXj)   =  ^laK\    »  (42) 

i=l         k=o 

where  f(x)  is  a  polynomial  of  degree  N-l,  to  be  chosen  in  some 
convenient  manner. 

We  want  to  determine  the  y.'s,  and  thus  we  require  that 
f(x)  =  f  (x)  such  that 

f.(x.)  =  6..  =f?  l  *  1    '  (43) 

an  orthogonality  condition.      If 

Vx)  -2CSq  **  •  <A4) 

J  K=0  J 

equation  (42)  reduces  to 

N-l 


y.  =  >   a  q  .     ,  (45) 

J    K=0       J 

where  q  „,  K  =  0,1,..., N-l,  are  determined  by  the  orthogonality 
condition  in  equation  (43) . 

Since  the  interval  of  integration  involved  in  the  preceding 
development  is  (0,1) ,  we  will  use  the  shifted  Legendre  polynomials 
in  approximating  the  integral  in  equation  (32).   Because  the  x.'s  are 
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the  zeros  of  the  shifted  Legendre  polynomial  P   (x) ,  we  can  use  the 

Lagrange  interpolation  formula 

PN*(x) 

The  desired  q._.'s  in  equation  (45)  are  the  coefficients  in 
this  polynomial  of  degree  N-l.   Repeating  this  procedure  for  each  j 
yields  the  desired  inverse  matrix  (q  .). 

INTERLACING  OF  THE  ROOTS 
It  should  be  noted  that  the  roots  (x.)  for  different  N  inter- 
lace (Appendix  B) .   Consequently  the  values  of  time 

t.  =  -log  x.  (47) 

l      °     l 

intertwine,  and  the  results  for  one  N  may  be  combined  with  those  of 
another  N  to  obtain  both  a  check  on  the  results  and  an  enlarged  set 
of  t  values  for  extrapolation  purposes. 


CHANGE  OF  TIME  SCALE 

The  method  of  inversion  described  for  the  Laplace  transform 
provides  approximate  values  for  the  function  u(t)  at  discrete  values 
of  t.   For  general  N,  we  obtain  u(t.)  at  the  points  t.  =  -log  xs  where 
x  are  the  roots  (or  zeros)  of  the  shifted  Legendre  polynomials,  V      (x) , 

It  is  obvious  that  the  inverse  transform  is  thus  determined  for 
a  very  small  range  of  t  values.   Since  increasing  N  does  not  signifi- 
cantly increase  the  range  of  t,  we  must  resort  to  another  method  of 
time  scaling.   Recall  that  if 

F(s)  =  L{u(t)}  = 


r°°  - 


e  Stu(t)dt,  (48) 


o 
then  from  Theorem  V,  it  follows  that 


L{u(at)}  =\   F(|j.  (49) 
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Therefore  using  the  approximate  relation 

£   A^'Va  log  x.)  -  \  F(f)  ,  (50) 

we  can  determine  u(t)  for  values  of  t.  =  -a  log  x„, 

ALTERNATE  METHODS 

Although  the  method  of  inverting  the  Laplace  transform  pre- 
sented involves  the  use  of  the  shifted  Legendre  polynomials,  there 
are  several  other  methods  that  can  be  used  in  obtaining  an  explicit 
inverse  matrix.  The  first  is  the  Method  of  Papoulis  which  utilizes 
the  Legendre  polynomials  and  the  second  is  based  on  a  different  set 
of  orthogonal  polynomials,  the  Laguerre  polynomials , 

Papoulis  Method 

The  Papoulis  method  of  inversion  involves  the  expansion  of  f(t) 
into  a  series  of  Legendre  polynomials  and  uses  the  classical  interval 
of  orthogonality  (-1,1).   It  is  clear  from  this  method  how  a  particular 
region  of  the  time  scale  may  be  examined  by  changing  the  scale  factor 

Y  [8]. 

Let  us  make  the  transformation 


x  =  e  Yt       y  >  0  (51) 


where  f(t)  becomes 


f(t)  =  f(-  ^  log  x)  E  f(x),  (52) 


Making  these  substitutions  into  equation  (1)  we  get 
YF(s)  = 


x(s/y)-l  £(x)dx  9  (53) 


J 
o 

By  expressing  s  in  terms  of  y  as  s  =  (2<+l)y,  equation  (53)  becomes 


YF(2ky+Y)  ■ 


f1    2k 


f(x)dx.  (54) 
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Defining  f(x)  as  an  even  function  (i.e.,  f(x)  =  f(-x)),  the  function 
f(x)  can  be  expanded  into  a  series  of  even  Legendre  polynomials  as 

00 

f(x)  =2  CNP2N(x)-  (55) 

N=o 

Using  the  fact  that  P„  (e  '  )  is  an  even  polynomial  in  e  ' 

— Yt 
of  degree  2N,  Papoulis  shows  that  the  transform  of  P?N(e   )  is  of 

the  form 

P(s)  =  P  (s)  =  (s-y)(s-3Y)...[s-(2N-l)T]  (56) 

nsj   *2lT8;     s(s+2y)   ...   (s+2Ny)  °b; 

Replacing  s  first  by  y   then  3y,  5y, . . . (2N+l)y  ...  in  equation  (56), 
gives  the  following  system  of  equations  to  evaluate  C  : 
yF(y)   =  Cq 

YF(3Y)  =  I  Co  +  ll  Cl 


(57) 


Ff2Nv+Y^  =  -H2_  +     2NCl      .      ,     2N(2N-2)  . . . 2CN 
jj^wy-nr;  -  2N+1  f  (2N+l)(2N+3)    "■    (2N+1)  (2N+3)  . .  .  (4N+1)  * 

Consequently  the  C  's  are  determined  and  the  function  f(x)  in  real 
time  is 

oo 

««)-2s  Ve"Yt>-  (58> 

N=o 

Gauss-Laguerre  Quadrature 
The  Laguerre  polynomials  are  the  unique  Gauss  quadrature 
polynomials  in  the  interval  (O,00)  with  weighting  function 

w(x)  =  e"ax.  (59) 

Denote  the  Laguerre  polynomials  by  L  (x)  where 

LN(x)  =  tt(x)(-1)N  ,  (60) 

and  the  interval  of  expansion  is  (0,°°),   These  polynomials  are  defined 

by  the  Rodrigue's  formula 

x  /  \    ok  d      /  N  -ax)  ,t,N 

Vax)  =  e   d^vx  e  /  '  (61) 
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It  follows  from  equation  (14)  that  the  orthogonality  condi- 
tion is 

•CO 

— rvv 
e    LN(ax)  LM(ax)dx  =  0       M  =h  N  (62) 


o 
Consequently  the  resursion  relations  are 

LN+1(x)  =  (1  +  2N  -  x)  LN(x)  -  N2LN_x(x)  ,  (63) 

and     LN+]/x)  "  (N+l)[L^(x)  -  LN(x)].  (64) 

The  least  squares  polynomial  approximation  to  f(x)  over  (O,00) 

with  respect  to  the  weighting  function 

w(x)  =  e"ax  ,  (65) 

N 
is      F(x)  =  X  bKVax>  »  (66> 

K=0 


where   b    =  -; — rr- 
k     (k!) 


00 

e"ax  f(x)  L  (ax)dx  .  (67) 


Consequently  the  coefficients  (or  weights)  become 

r°°  L„(ax) 

A.    =   e    7 sT  77 r  dx  .  (67) 

i     J      (x-x  )L^(axi) 

0 

Using  the  properties  stated  above,  a  similar  procedure  to  that 
employed  with  the  Legendre  polynomials  may  be  performed  and  an  explicit 
inversion  formula  derived. 

TEST  CASES 

A  numerical  technique  has  been  developed  based  on  shifted 
Legendre  polynomials  to  invert  the  Laplace  transform  of  a  function. 
Let  us  now  examine  the  accuracy  obtained  in  utilizing  this  technique. 

Two  test  functions  are  considered.   The  first  function  is 
u(t)  =  constant  =  1  where  the  Laplace  transform  is  v(s)  =  —  .   The 

8 

second  test  function  is  u(t)  =  e   where  the  Laplace  transform  is 

v(s)  =  — pr  •   The  results  obtained  for  these  test  cases  are  tabulated 
s+l 
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in  Tables  I  and  II.  Notice  that  the  accuracy  decreases  as  N  (number 
of  roots)  increases.   However  8-9  figure  accuracy  is  obtained  in 
test  cases  I  and  6-9  figure  accuracy  is  obtained  in  test  case  II. 
These  results  indicate  that  the  numerical  inversion  technique  used 
is  very  accurate. 
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TABLE  I 


TjStfT   CSB   I 


N 
4 


10 


14 


CX^CP   VALUE 

0.999999999D  00 
0.999999999D  00 
0.999999999D  00 
0.999999999D  00 

0.999999999D  CO 

C999999999D  00 

0.999999999D  00 

0.999999999D  00 

0.99.999999D  00 

0.999999999D  00 

0^-9999999991)  00 

0.99999999913  00 
0.999999999D  00 
0.999999999D  00 
0.999999999D  00 
0.999999999D  00 
0.999999999D  00 
J.999999999D  00 
0.999999999D  00 
0.999999999D  00 
0.999999999D  00 

0.999999999D  00 
0.999999999D  jO 
0.99^9999991)  00 
0.999999999D  00 
0.999999999D  CO 
0.999999999D  00 
0.999999999D  00 
0.999999999D  jO 
0.999999999D  00 
0.999999999D  JO 
0.v99999999D  00 
0.999999999D  00 
0.999999999D  00 
0.999999999D  00 


APPROXIMATE  VALUE 

0.100000000D  01 
0.999999999D  00 
0.999999999D  00 
0.999999999D  00 

0.999999999D  00 
0.100000000D  01 
0.100000000D  01 
0.9999999992  00 
0.100000000D  01 
0.999999999D  00 
0.100000000D  01 

0.100000000D  01 
0.999999999D  00 
0.100000000D  01 
0.100000000D  01 
0.999999999D  00 
O.lOOwOOOOOD  01 
0.999999999D  00 
0.100000000D  01 
0.999999999D  00 
0.100000000L  01 

0.999999966D  00 

0.999999962D  00 
0.100000012D  01 
0.100000003D  01 
0.999999988D  00 
0.1000000C1D  01 
0.999999891D  00 
0.999999944D  00 
0.100c00004D  01 
0.999999862D  00 
0.100000009D  01 
0.999999977D  00 
0.999999962L  00 
0.100000005D  01 
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TABLE  II 

PEST    CAoi^   II 


N 
4 


7 


10 


14 


EXACT    VALUE 

O.930568155D  00 
0.66999052113  00 
0.330009478D  00 
0.694318442L-01 

u. 97455 3956D  00 
O.o707b559213  00 
0.70292257513  00 
0.499999999D  00 
0.^9707742413  00 
0.129234407T  00 
0.254460438D-01 

u.96o953264D  00 
0.9325316633  00 
0.0397047641  00 
0.71669769713  00 
0.^7443716913  00 
0.425;b2d30D  00 
0.2633023021)  00 
0.160^9521513  00 
0.674663166D-01 
0.1304o7357D-01 

0.99314190413  00 
0.  .6421744113  00 
0.9136^065713  00 
0.843646452D  00 
0.757624318D  00 
0.65955bl84D  00 
0.554027474'D  -0 
0.445972525D  00 
0.340443815D  00 
0.242375681D  00 
0.15635354713  00 
0.663993424D-01 
O.357o2558ll3-Ol 
0.6b5b09565D-02 


APKiOXIMA  i?£    '/ALUE 

0.930568155D  00 
0.669990521D  00 
0.330009478D  00 
0.694318442D-01 

0.974553956D  00 
0.870765592D  00 
0.702922575D  00 
0.500000000D  00 
0.29707742413  00 
0.129234407D  00 
0.254460438D-01 

0.986953264D  00 
0.932531683E  00 

0.83970478413  00 
0.71669769713  00 
0.574437169D  00 
0.425562830D  00 
0.283302302D  00 
0.160295215D  00 
0.o74b631b3L-01 
0.130467358D-01 

0.99314183313  00 
0.964217450D  00 
0.913600761D  00 
0.643646317D  00 
0.757624350D  00 
0.65955616513  00 
0.55402738913  JO 
0.44597258713  00 
0.3404^379313  00 
0.242375649D  00 
0.156353551T3  00 
0.863993391D-01 
0.357825607D-01 
0.685815140D-02 
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III.   THE  INVERSE  PROBLEM  IN  TRANSIENT  HEAT  CONDUCTION 

INTRODUCTION 
The  problems  presented  in  this  section  are  representative  of 
typical  one,  two,  and  three  dimensional  geometries.   They  were  selected 
because  of  the  existence  of  exact  solutions,  necessary  for  comparison 
of  results.   However,  more  complicated  geometries  can  be  accommodated 
by  this  method  of  solution,  but  the  equations  and  analysis  become 
quite  complex.   The  FORTRAN  IV  programs  used  in  the  solution  of  these 
problems  are  listed  in  Appendix  A. 

DIMENSIONLESS  VARIABLES 
The  solution  to  a  problem  in  heat  conduction  can  always  be 
expressed  in  terms  of  dimensionless  variables.   This  is  usually  done 
when  making  numerical  calculations  and  is  adopted  here  because  of 
simplification  in  the  differential  equations  involved  and  generality 
in  application  of  the  solutions  obtained. 

The  dimensionless  variables  used  in  the  following  problems  are 

Kt   Kt 

T     £2    R2 

where   K  = 


Cp (density)  ' 

x   _r 

I   ~'~~   R 


uK    uK 
u  = 


6£2   3R2 


%  =  R  YN 
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THE  SEMI-INFINITE  SOLID 

The  method  of  numerical  inversion  of  the  Laplace  transform 
is  now  applied  to  the  solution  of  the  semi-infinite  solid.   Since 
there  exists  no  physical  characteristic  length,  dimensional  variables 
are  used  in  the  solution  to  this  problem. 

The  partial  differential  equation  governing  the  temperature 
function  is 

with  the  boundary  conditions 

u(0,t)  =  Bt  ,  (70) 

and 

u(x,0)  =  0   .  (71) 

The  solution  of  this  direct  problem  is 

u(x,t)  =  4Bti2erfc  ~—     .  (72) 

J2  Kt 

Using  this  solution  as  the  known  temperature  function  at  x.  (an 

interior  location),  we  can  solve  the  inverse  problem  to  obtain  u(0,t) 

at  x  =  0  (the  surface  of  the  solid).   Here  the  exact  solution  is  known 

as 

u(0,t)  =  Bt  .  (73) 

Transforming  the  differential  equation  into  the  s  plane  we  obtain  the 

ordinary  differential  equation 

0-fv-O  (74) 

where  v  is  the  Laplace  transform  of  u(x,t)  and  is  a  function  of  x  and 
s.   The  solution  to  equation  (74)  is  of  the  form 

v(x,s)  =  A  e"^  X  (75) 

Solving  this  equation  at  x  -   x.  we  find  that 

A  -  v(x1,s)e>'s/K  Xi  (76) 
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and  therefore  the  solution  of  equation  (74)  becomes 

/   \     /    .  ^- Js/K  (x-Xt )  ,.,-.. 

v(x,s)  =  v(x-,s)e  N       1  (77) 

Since  the  temperature  function  at  an  interior  location  x1  is  given 

in  equation  (72),  we  can  solve  equation  (77)  by  taking  the  Laplace 

transform  of  equation  (72)  which  then  becomes 


v(xn  ,s)  =  — 5-  e  M     1  (78) 


Substituting  equation  (78)  into  equation  (77)  we  have  the  solution 

in  the  s  plane  as 

,  .        /B        -JsTk  x.)    -IsTk   (x-x.)  ,,on 

v(x1,s)    -  f-y  e  >l  ye  si  V  (79) 

At   the  surface  of   the  solid  x  =   0,    and   the  solution  is 

v(0,s)   -  K  e'sF7^  xl)eJ^  xl  (80) 

Taking  the  Laplace  inverse  of  equation  (80)  we  obtain  the  solution 

in  the  time  plane  at  x  =  0  (the  surface  of  the  solid). 

An  illustration  of  the  semi-infinite  problem  is  presented 

in  Figure  1.   In  obtaining  the  results  plotted  in  Figure  2  it  was 

found  advantageous  to  use  a  convergence  factor.   The  convergence 

factor  e   was  added  to  the  function  u(x,t)  to  improve  the  accuracy 

of  solution.   The  function  u(x,t)  was  multiplied  by  e   ,  transformed 

into  the  s  plane  and  subsequently  inverted  back  into  the  time  plane. 

The  final  answer  was  then  multiplied  by  e  .  Note  that  the  Laplace 

transform  of  this  function  is 

L{e_tu(x,t)}  =  F(s+1). 

THE  LONG  CYLINDER 
The  basic  assumption  involved  in  the  solution  of  the  "long" 
cylinder  is  that  the  cylinder  is  sufficiently  long  so  that  end  effects 
can  be  neglected.   In  other  words  conduction  is  considered  to  be  only 
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in  the  radial  direction  (Figure  1) .   The  governing  partial  differential 
equation  for  the  temperature  function  u(r,t)  is 

subject  to  the  boundary  conditions 

u(R,t)  =  St  ,  (82) 

and     u(r,0)  =  0   .  (83) 

In  terms  of  dimensionless  variables,  the  solution  of  the  direct 

problem  can  be  represented  as  the  sum  of  an  infinite  series  as 

where  the  a^   are  defined  by 

JQ(aN)  =  0   .  (85) 

Using  this  solution  as  the  known  temperature  function  at  p. 

(an  interior  location) ,  we  can  solve  the  inverse  problem  to  obtain 

u(p,x)  at  p  =  1  (the  surface).  Here  the  exact  solution  is  known  as 

u(1,t)  =  T   .  (86) 
Proceeding  as  before  we  find  that  the  solution  of  the  inverse 

problem  in  the  s  plane  is 

v(l,s)  =  v(p,,s)  S"^-'  \  .  (87) 


^rs)i^f^   > 

T(p  ,.)  ■  ^  -  J^ll  +  2  ±  U-\       l°W*>    .  (88) 

KiJV  s2    4s      §=I^N2+s/  otN3  Jx(aN) 


where  ,    i  ^  2 


By  numerically  inverting  equation  (87) ,  we  obtain  the  solution  of 
the  inverse  problem  at  the  surface  (p  =  1)  where  the  exact  solution 
is  given  by  equation  (86) . 

The  results  obtained  for  the  "long"  cylinder  are  plotted  in 
Figures  5  and  6.   It  should  be  noted  that  as  the  distance  from  the 
surface  is  increased,  the  accuracy  of  solution  decreases  (Figure  5). 
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This  is  reasonable  since  the  temperature  will  be  dampened  more  as  the 

distance  is  increased.   Let  us  consider  a  characteristic  parameter 

1-r 

2jKt 
If  a  thermal  penetration  depth  is  taken  as 

6  =  2jlCt   , 

then  physical  significance  can  be  given  to  r\.      Now  we  can  examine  the 

accuracy  of  solution  as  r\   is  varied.   In  Figure  5,  surface  temperature 

is  plotted  versus  n  for  6  =  constant.   It  can  be  seen  that  as  r\   is 

increased,  the  accuracy  of  solution  decreases.   Since  the  only  variable 

in  n  is  now  r,  it  is  obvious  that  as  r  increases,  the  accuracy  of 

solution  increases.   By  interpreting  small  r\   as  "close"  to  the  surface 

and  large  n  as  "far  away"  from  the  surface,  we  can  see  the  effects 

of  damping  on  the  solution. 

A  convergence  factor  was  also  used  in  this  problem  to  improve 

the  accuracy.   As  before,  e   was  used  as  the  convergence  factor,  and 

it  can  be  seen  from  Figure  4  that  the  solution  is  improved.  However 

the  results  obtained  are  not  quite  as  good  as  the  other  problems 

considered.   The  inaccuracy  introduced  into  this  problem  is  due  to 

the  evaluation  of  an  infinite  series  of  Bessel  functions. 

THE  SPHERE 
The  basic  assumption  involved  in  the  solution  of  the  sphere 
is  axial  symmetry.   In  other  words  the  temperature  function  is  not 
dependent  on  $  (Figure  1) .   The  governing  partial  differential  equa- 
tion for  the  temperature  function  u(r,t)  is 


?-*(lf+f£).  <*« 
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subject   to   the  boundary  conditions 

u(R,t)    =   Bt      ,  (90) 

and 

u(r,0)    =  0        .  (91) 

In  terms  of  dimensionless  variables,  the  solution  of  the 
direct  problem  can  be  represented  as  the  sum  of  an  infinite  series 
as  oo      n 

"  r  N         1-p2     2   ^T"  (-1)    -N27T2T  .  ...       N       .... 

u  (p,T)  =  T  -  -f-  -  -T-  2L   "Sp-  e      sm(NTrp).     (92) 

M  N=l 
Using  this  solution  as  the  known  temperature  function  at 

p   (an  interior  location) ,  we  can  solve  the  inverse  problem  to 

obtain  u(p,T)  at  p  =  1  (the  surface).   Here  the  exact  solution  is 

known  as 

u(l,x)  =  T  .  (93) 

Proceeding  as  before  we  find  that  the  solution  of  the  inverse 

problem  in  the  s  plane  is 

P]_  sinh  (sfs) 
v(l,s)   =  v(Pl,s)  sinh  q=  pp   .  (94) 

where  ±  1_pi2    2 


oo  N  /  \ 


v(p,  ,s)  =  —5-  - 

1      s     6s     ,n"3p1  ^ 

By  numerically  inverting  equation  (94)  ,  we  obtain  the  solution  of 

the  inverse  problem  at  the  surface  (p  =  1)  where  the  exact  solution 

is  given  by  equation  (93) . 

The  results  obtained  for  the  sphere  are  plotted  in  Figure  7. 

Here  as  before  a  convergence  factor  e   was  used  to  improve  the 

solution.   The  results  shown  in  Figure  7  are  in  very  good  agreement 

with  the  known  exact  solution. 
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INVERSE    PROBLEMS  IN  TRANSIENT 
HEAT  CONDUCTION 
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SURFACE  TEMPERATURE  vs.  TIME 
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SURFACE  TEMPERATURE  vs.  TIME 
FOR 
"LONG"  CYLINDER 
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DISCUSSION  OF  RESULTS 

In  the  determination  of  surface  temperature  by  the  method  of 
numerical  inversion  of  the  Laplace  transform,  it  was  found  necessary 
to  utilize  a  convergence  factor  to  improve  accuracy.   The  function 
u(x,t)  was  multiplied  by  e   ,  transformed  into  the  s  plane,  and 
numerically  inverted  into  the  time  plane.   Then  the  function  obtained 
was  multiplied  by  e  to  obtain  the  desired  surface  temperature  function 
u(0,t).   This  method  caused  the  approximation  to  converge  to  the  known 
exact  solution.   In  the  semi-infinite  solid  and  the  sphere,  the 
results  obtained  from  the  use  of  the  convergence  factor  were  very 
good.   In  the  "long"  cylinder  the  results  obtained  were  not  quite  as 
good.   The  inaccuracy  introduced  in  this  problem  is  due  to  the  evalua- 
tion of  an  infinite  series  of  Bessel  functions. 

The  damping  effects  were  considered  for  the  case  of  the  "long" 
cylinder.   It  was  found  that  as  the  distance  from  the  surface  was 
increased,  the  accuracy  of  the  solution  decreased.   This  was  a 
reasonable  result,  since  it  was  logical  to  expect  the  solution  to 
lag  as  the  distance  from  the  surface  was  increased. 

The  use  of  time  scaling  (mentioned  in  section  II)  was  only 
incorporated  in  the  semi-infinite  solid  FORTRAN  IV  program,  but  it 
can  be  included  in  the  other  problems .   Although  the  initial  tempera- 
ture distribution  was  taken  as  zero,  any  reasonable  initial  tempera- 
ture distribution  can  be  handled  with  this  method .   Such  a  change 
would  simply  result  in  more  terms  to  be  inverted.   Also  it  should  be 
noted  that  although  the  solution  was  obtained  for  the  surface  tempera- 
ture, heat  flux  at  the  surface  can  also  be  determined. 
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CONCLUSIONS 

Three  geometries  were  considered  in  the  solution  to  the 
inverse  problem  in  transient  heat  conduction.   These  geometries, 
the  semi- infinite  solid,  "long"  cylinder,  and  sphere,  are  typical 
of  one,  two,  and  three  dimensional  problems  respectively  (it  should 
be  noted  here  that  an  analogous  approach  can  be  made  in  the  solution 
to  the  finite  plane  bounded  by  two  parallel  planes) .   Solutions  were 
obtained  which  agreed  very  well  to  the  known  exact  solution. 

Several  parameters  were  investigated  with  respect  to  the 
effects  on  the  accuracy  of  solution.   It  was  found  in  all  cases  that 
utilizing  a  convergence  factor  caused  the  solution  to  better  approxi- 
mate the  known  exact  solution.   The  damping  effect  was  considered  in 
the  "long"  cylinder  problem.   As  expected  it  was  found  that  as  p, 
increased  (distance  from  the  surface  decreased)  the  accuracy  of 
solution  increased. 

Therefore  the  solutions  obtained  indicate  that  the  method  of 
numerically  inverting  the  Laplace  transform  provides  an  accurate 
approximation  to  the  known  exact  solutions. 


42 


IV.    SUMMARY 

The  numerical  technique  involved  in  numerically  inverting  the 
Laplace  transform  was  methodically  developed.   Elementary  theorems 
of  the  Laplace  transform  were  introduced,  and  the  instability  of  the 
inverse  of  the  Laplace  transform  was  briefly  discussed.   Gaussian 
quadrature  methods  were  subsequently  introduced  and  developed.   The 
shifted  Legendre  polynomials  were  then  examined  for  several  reasons. 
It  was  found  that  the  coefficients  of  the  shifted  polynomials  were 
integers,  and  the  region  of  orthogonality  (0,1)  was  found  to  be  most 
convenient  in  the  following  development.   Then,  an  explicit  inverse 
matrix  was  found  which  would  be  used  later  to  solve  the  inverse  prob- 
lem in  transient  heat  conduction 0   Two  useful  numerical  properties 
were  discussed.   Interlacing  of  the  roots  (and  therefore  times)  was 
found  useful  in  interpolation  of  the  results,  while  change  of  the 
time  scale  was  found  useful  in  extending  the  range  of  times.   Two 
alternate  methods  of  solution  (Papoulis  and  Gauss-Laguerre  quadrature) 
that  could  be  used  to  obtain  an  explicit  inverse  matrix  were  subse- 
quently introduced.   Finally,  the  explicit  inverse  matrix  was  tested 
for  accuracy.   The  two  test  cases  employed  were  for  the  functions 
u(t)  =  1,  and  u(t)  =  e   .   The  accuracy  was  found  to  be  very  good  for 
both  cases. 

Three  geometries  were  considered  in  the  solution  of  the 
inverse  problem  in  transient  heat  conduction.   These  problems,  the 
semi- infinite  solid,  "long"  cylinder  and  sphere,  are  typical  of  one, 
two,  and  three  dimensional  problems  -   They  were  chosen  because  of  the 
existence  of  exact  solutions.   It  should  be  noted  that  the  results 
obtained  were  for  discrete  points  and  not  for  a  continuous  spectra  of 
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information.   Consequently  this  method  is  applicable  to  problems 
where  discrete  information  is  known  (i.e.,  data  from  experiment). 

Several  parameters  were  examined  with  respect  to  influence 
on  the  accuracy  of  solution.   A  convergence  factor  e   was  found 
to  improve  the  accuracy  of  solution  in  all  cases  considered.   It 
was  found  that  as  the  distance  from  the  surface  was  increased, 
the  accuracy  of  solution  decreased.   This  was  considered  reasonable 
since  damping  effects  should  be  more  noticeable  at  greater  distances 
from  the  surface. 

Consequently  the  solution  to  the  inverse  problem  in  transient 
heat  conduction  by  numerical  inversion  of  the  Laplace  transform  was 
found  to  yield  accurate  results  for  semi-infinite  solid,  "long" 
cylinder,  and  sphere. 
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APPENDIX  A 
FORTRAN  IV  PROGRAMS 

1.   The  Semi-infinite  Solid 
The  following  program  is  listed: 

SEMI-INFINITE  PLANE  SOLUTION  USING  DIMENSIONLESS 
VARIABLES  IN  ONE-DIMENSIONAL  SYSTEM, 

The  following  subroutines  are  needed: 
INPUT  Subroutine 

EXPINV  Subroutine 
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FORTRAN  IV  PROGRAM  SYMBOLS 
MAIN  PROGRAMS 

A  Explicit  inverse  matrix 

R  Roots  of  the  shifted  Legendre  polynomials  (x„) 

W  Weights  of  the  shifted  Legendre  polynomials  (A„) 

Tl  Values  of  time  (t,  =  -log  x.) 

Kt   Kt 

TAU  Non-dimensional  time  (t  =  -rj  =  ^y) 

RO  Distance  (p) 

R01  Interior  location  (p..) 

APPR  Approximate  value  of  v(x,s) 

XK  Reciprocal  of  K 

EXACT  Known  exact  solution 

SCALE  Scale  factor  (a) 

C  Value  of  e~  u(p,x) 

ALFA  Roots  of  the  equation  J  (cO  =  0 

ARG  Argument  of  the  function  J  (aNp-^) 

ARQ  Argument  of  the  function  I  (Js  p..) 

ARP  Argument  of  the  function  I  ( fs   ) 

BESSO  Resultant  value  of  J  (o^P-^ 

BESS1  Resultant  value  of  ^-,(0 

BISS  Resultant  value  of  I  (JT  ) 

o  ^ 

BISS1  Resultant  value  of  I  (Js  p.,) 

ADA  Value  of  n, 

N  Number  of  roots 
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EXP1NV  SUBROUTINE 

C      Legendre  coefficients 

B      Coefficients  of  the  polynomials  P  (x)/(x-x„) 

P      Legendre  derivatives 
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c 
c 
c 
c 
c 
c 
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SOLUTION  OF  THE  INVERSE  PROBLEM  IN  TRANSIENT  HEAT  CONDUCTION 

BOUNDARY  CONDITIONS  ARE: 

. .  ...   1.  INITIAL  TEMPERATURE  OISTRIRUTTON  IS  ZERO. 

2.  SURFACE  TEMPERATURE  DISTRIBUTION  IS:   U(1,TAU)*TAU 

I.  SEMI-INFINITE  PLANF  SOLUTION  USING  DI MENSIONLESS  VARIABLES 

—  iN-ONE-DIMENSIONAL SYSTEH, 


*R 


*R 


20 


24 


*E) 


DIM 

( 15), E 

- — REA 

(20)  ,R 

00 

CAL 

CAL 

WRI 

DO 

WRI 

SCA 

XK  = 
PI  = 
00 

TAU 

IP  = 
R01 
IR  = 

R0( 

CON 

DO 

Q=I 

S=0 

EPP 
APP 


IN  T1(15).A(15,15),R(15),W(15),TAU(15),C(15),RQ(11)?APP 

15).F(20) 

•» 

) 

U=2,15 
'UT(N,R,W,TD 

INV(N-r^iWtA) 

106) 

1,M 

102)R( I ) ,W( I ) 

CD*00 

C2 

59265 

1«N 

<C.OD«-00 

,00*00 

i=l,N 

F 

=DEXP(-DSQRT(S*XK)*R01( IP) )/S**2 

=EPPR( IS)*DEXP(-DSQRT(S*XK)*(RO( IR)-R01< IP) ) )/(XK*SCAL 


ENSIOI 

XACTI 

L*8  A,RTWtC^Tl^T^^r4APRtEXACT,F,XKvR0TAT)A(^0VyS^SCAL^rXH€-P-P— 

01  (20 

100  J 

L  INPI 

L  EXP 

TE(6, 

20  1  = 

TE(6, 

LE=2. 

l.OD  +  i 

3.141 

4  I  1  = 

(  I  U  =  SCALE*T1(  H)/XK 

1 

(  I  P  I 

1 

IR)=0- 

TINUF 

40  IS  = 

S+l 

/SGAL 

R(IS) 

R(  IS) 


40 


61 
-62- 


400- 


-EXA 
ADA 
CON 
DO 

-C(  I 
DO 
C(I 
DO 

-F(J 
WRI 
WRI 
WRI 

-WRI 
IP  = 
R01 
IF( 

-CGN 


N 

■00- 
N 


N 


60 

j 

—70- 

80 

83 

102 

106- 

108 


CAL 
FOR 
20  X, *A 

FOR 

FOR 
FOR 
FOR 

FOR 

FOR 
END 


CT(  IS)=TAU( IS) 

( IS)=P01( IP)/(2.0D*00*DSQRT(T1(IS)/XK)) 

TINUE 

61  1=1, 

)=O.CC+ 

61  J=li 
)  =  C(  I)*A(  I,J)*APPR< J) 

62  J=l, 
)=C(J)* 
TE(6,83 
TE(6,80 
TE<6,60 

TE  (  6  ,  70 )  (T4-m^M4-K£X  AC 
IP+l 
( IP)=RO 

roh  ip) 

TINUE 

L  EXIT 

MAT(//,11X, "TIME*, 14X, 'APPROXIMATE  VALUE ', 15X ,' EXACT  VALUE* 

DA» ♦//) 

MAT(4(2Xt-D2A.L7-)»/)- 

MAT(//, 

MAT(//, 

MAT(//, 

MAT(//r 

MAT(//, 


iDEX-P-l-J-UJ-W- 

)N 

)R01( IP) 

) 


1(  IP-1H0. 10*00 
.LT. 1.000)  GO  TO  24 


H~lrrW- 


3X,«THE  VALUE  OF  R01  IS   «,F6.4,//) 
3Xf'THE  VALUE  OF  N  IS   ••15,//) 
2(2X,D24.17) ,/) 

10X^ 'ROOTS  «^rl3Xt  'WEIGHTS'^//) 

5X,'T  a-LOG  R«,20X,/f 024.17,/) 
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2.   The  "Long"  Cylinder 
The  following  program  is  listed: 

TWO-DIMENSIONAL  "LONG"  CYLINDER  SOLUTION. 

The  following  subroutines  are  needed: 

INPUT  subroutine 

EXPINV  subroutine 

BESI  subroutine 

BESJ  subroutine 
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SOLUTIHN  OF  THE  INVFRSE  PROBLEM  IN  TRANSIENT  HEAT  CONDUCTION 

BOUNDARY  CONDITIONS  ARE: 

1,  INITIAL  TEMP  DISTRIBUTION  IS  ZERO. 

2.  SURFACE  TFMP  DISTRIBUTION  IS  :  U  (  1  ,  T  AU  )  =T  All. 
III.  TWO-DIMENSIONAL  "LONG"C YL I NDER  SOLUTION. 

DIMFNSION  TKlS)fA(15,15),R!15),W!l5),TAU!l5),C!15),Rn!15),APP 
*R(15).FXACT(15),ARG(21),SUMI15KALFA(15)fARP(2O>,ARO(2r>),ROl(?0» 

PFAL*R  A,R,W,C,T1  ,TAU,APPR,EXACT, ALFA, SUM, BESSO, BESS1 , B ISS. B I S 
*SltARG,ARQ,ARP,S,XK,F<201 

no  2  i=l,in        

2      RFAD(5,1 in)ALFA( I ) 

DO  100  JJ=2,15 

CALL  INPUT(N,R,W,T1) 

CALL  EXPINV(N,R,W,A) 

WRlTE(6t 1C6) 

DC  20  I=1,N 
20      WRITE!6,102)R! I),W(I) 

XK=1. 60*01  -   

DO  4  I  1  =  1, N 
4      TAU! II)=T1! I I)/XK 

PI  =  3. 14159265        ._ 

R0?IR)=1. 00*00 

IP=1 

ROK  IP)=0.C0*00 
24      CONTINUE 

DO  90  IS=1,N 

SUM< IS)=0.0D*0C 

S  =  I  S ♦! 

DO  6  K=1,1C  •  —         -  

ARG(K)=R01( IP)*ALFA(K> 

CALL  BESJ(ARG,0,BESS0,1.0D-12,IER) 

CALL    BESJ(ALFA,1,BESS1.1.0D-12,IER) 
6  SUM( IS)=SUM! IS)*!  1.00*00/ < ALFA! K )**2*S ) >*BESSO/< ALFA! K )**3*BES 

*S1) 

ARP(  IS)=DSORT( S) 

ARO(IS)=R01( IP)*DSQRT!S) 

CALL  BESI(AR0,C,BISS1,IER>   - - 

CALL  BESI(ARP,0,BISS. IER) 

APPRUS)=(  t  1.0D+00/S**2)-!<1.0D*00-R01(IP)**2)/<4.0D«-00*S))*2. 
*OD+00*SUM( IS) )*BISS/(BISS1*XK) 

EXACT! IS)=TAU( IS) 
90      CONTINUE 

DO  50  1=1, N 

Cf I )=O.OD+00 

"  .  DO  50  J=1,N  --  

5C      C(I )=C< I)*A( I,J)*APPR!J) 

DO  62  J=1,N 
62      F(J)=C! J)*PEXP!T1!J) ) 

WRITE(6,83»N 

WRITE(6,B0)R01( IP) 

WRITE! 6,6C  ) 

WRITE  (6,  70)  (TAU!  I  ),F(  I),  EXACT  (  I),  I*  UN) 

ip=ip+1  ......  

ROK  IP)=R01(  IP-1)*0.1D*00 
IF(R01( IP). LT. 1.000)  GO  TO  24 
100      CONTINUE 
CALL  EXIT 
60      FORMAT!//, UX, 'TIME' ,14X, 'APPROXIMATE  VALUE* , 15X, • EXACT  VALUE* 

*,//) 
70      F0RMAT(3(2X,D24.17)  ,/) 

80  -—  FORMAT!//, 3X, 'THE  VALUE  OF  R01  IS   •  ,F6.4.  //  >      ~ 

83      FORMAT!//, 3X, 'THE  VALUE  OF  N  IS   «,I5,//) 
1C2      FORMAT!//, 2!2X, 024. 17) ,/) 
106      FORMAT!//, 10 X, •ROOTS* ♦ 13X, • WE IGHTS* ,  // ) 
1C8      FORMAT!//, 5X,«T  «-LOG  R  •  ,20X, /,D24. 17, / ) 
110      F0RMAT!013.6) 
ENO 
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3.   The  Sphere 
The  following  program  is  listed: 

THREE-DIMENSIONAL  SPHERE  SOLUTION. 

The  following  subroutines  are  needed: 
INPUT  subroutine 

EXPINV  subroutine 
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SOLUTION  OF  THE  INVERSE  PROBLEM  IN  TRANSIENT  HEAT  CONDUCTION 

BOUNDARY  CONDITIONS  ARE: 

1.  INITIAL  TEMPERATURE  DISTRIBUTION  TS  ZERO. 

2.  SURFACE  TEMPERATURE  DISTRIBUTION  IS:   U(1.TAU)=TAU 

IV.  THREE-DIMENSIONAL  SPHERE  SOLUTION--- -■   


DIMENSION  T1(15),A(1«>,15),R<15).W(15) »TAU<15) »C( IS) f RO( 1 1 ) • APP 

*R (15) , EXACT!  15), SUM!  15) , ROl! 21) 

SEAL* 8  A,R,W,C,T1,  TAU, APPR ,E XAC T, SUM, F ! 20 ) 

DO  100  JJ  =  2,  15 

CALL  INPUT(N,R,W,T1) 

CAl L  EXPINV(N,R,W,A) 

_   WRITE (A*  106)       

DO  20  1=1, N 
20      WRITF(6,  102)R(  I) tW( I ) 
XK= 16 . OOCOOOOOOOOOOOO 
DO  4  I  1  =  1, N 
4      TAU( I  I  )  =  T1(  I  I  ) /XK 
PI=3. 14159265 
IR=1 

ROI IR)=1. 00+00  

IP=1 

R01I IP)=O.OD+00 
24      CONTINUE 

On  go  IS=i,N 
SUM( IS)=C.0D+OO 
00  6  K=l,7 
S=IS+1 

— 6 SUM(  I  S)=-SUM*  l_S)+r(-i)**K/K**3)*-(  K00+00/!K**2*PI**2  +  S)  )*SIN<K*- 

*PI*RQ1(  IP)  ) 

APPR(  IS)  =  (  (  1.0D  +  00/S**2)-<  (  1.0D+00-R01!  I  P  )  **2)  /  (  6.0D«-00*S  )  )-(2 
♦.00+00/ ( PI** 3*R01( IP) ) )*SUM( IS) ) ♦ROl ( IP) *SI NH( SORT( S ) ) /< SI NH( SORT! 
*S)*R01  (  IP  )  )*XK)  ...   - 

EXACT( IS)=TAU( IS) 
90      CONTINUE 

DO  50  I  =  1,N 

C( I )  =  0.00 +00 

DO  50  J=1,N 
50      C( I >=C(  I)+A(  I, J)*APPR( J) 

DO  62  J=1,N 
62      F(J )=C( J )*DEXP(T1( J) ) 
WRITE!6,8C)R01( IP) 
WRITE(6,60) 

WRITE (6, 70) (TAU! I ),F! I ), EXACT! I ),I=1,N) 
jp^IP+1 _ _ 

ROl  (  IP  )=P01(  IP-D+0.1D+00 
IF(R01( IP) .LT. 1.000)  GO  TO  24 
100      CONTINUE 
CALL  EXIT 
60      FORMAT!//, 11X, 'TIME' ,14X, 'APPROXIMATE  VALUE •» 15X ,' EXACT  VALUE1 

♦  ,//) 
70      F0RMAT(3(2X,D24.17) ,/) 

80 F0RMAT(//^3X,  'THE  VALUE- -OF— R-ai—l-S *^F  6»  4  ,  /  /  ) 

102      FORMAT!//, 2!2X,D24. 17) ,/) 

106      FORMAT!//, 10 X, 'ROOTS' ,13X, 'WEIGHTS' ,//) 
108      FORMAT!//, 5Xf «T  =-LOG  R • ,20X, /, 024. 17,/ ) 
END — —  —   
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4.   Subroutines 
The  following  subroutines  are  listed: 
INPUT  subroutine 

EXPINV  subroutine 

BESI  subroutine 

BESJ  subroutine 


* 

"Subroutine  BESJ"  and  "subroutine  BESI,"  System/360 

Scientific  Subroutine  Package  (360A-CM-Q3X)  Version  II  Programmer's 
Manual.   White  Plains,  New  York:   International  Business  Machines 
Corporation,  1967,  pp.  156,  157,  159. 
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10 

20 

15 

40 

SO 
50 


Sl.'llRniJTINF     FXPINV(N,P,W,A) 

HI  VFNSIflN     P(15l, W(lSV,A(]f>, IS),  f.  (17,17),  P(l^',  lSI.PCim 

RfAL*H    R.W,A,C»Xl,X2,X,XI,R,P7FRn,PKFV?,PP(Vl,/l,7?,73,PPFV3,P 

*,ni:\r.v 

LECENOPF  COEFFICIENTS 

IS! PI  |)S  =  N+  1 

no  10  m=i.npu»s 
C  (  m  ,  1 )  =  1  .  o  i  -  ♦■  n  "j 
on  ?o  k=i.nplus 

OG  20  M=lf|MPLUS 

IF(K.f,T.M)   C( M,K) =0.0  0+00 

CCNTIMIIF 

00  40  K=1,NP|_US 

nn  ao  v=i,mphjs 

IF(K.GT.M)   on  TO  40 


*y. 


C(M  +  ltK  +  l)  =  ((2*M-l)*C(MfK4-l)-(M-n*C(M-l,K*l 

CONTINUE 


)-(4*M-?)*C(M,K ) )/ 


B  COEFFICIENTS 

•-0  0  8C  I=1,N 

P( I ,N)=C(N+1,N+1) 

CONTINUE 

00  90  1=1, N 

XI  =  RU  ) 

K  =  N 

K  =  K  —  1 

B(I,K)=C(N+1,K+1)+XI*B(I,K+1) 

IF(K.NE.  1)  GO"  T0--50^ 

CONTINUE 

LEGENORE  OFRIVATIVES 


150 
180 


NWINUS=N-1 
00  180  1=1, N 
X  =  R  (  I  ) 

P7FRO-l.OO--0  0-2.0n-<7C«*X 

PRFV2=P7FR0 

PPEV1=1.00  00 

OC  150  M=2,NMINUS 

7  2=2*M-1 

73=V 

PRFV3=-( 71*PREV1-Z2*P7ER0*PREV2)/Z3 

PRFV1  =  PPEV? . ; 

PRFV2=PRFV3 
Z1  =  N 


7  2=2.on  or-*x*(x-i.cn  oo) 

Pi  I  )=PREV3*Z1/Z2 
EXPLICIT  INVERSE  MATRIX 


22C 


-nn-220  -m=i-,n 

OFNn-V  =  W(  V)*0  (M) 
00    220    K=1,N  v 
A(M,K)=B(M,K  )/vOENOM 
RETURN 

END 
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10 

11 

12 
100 

1 G  1 


S1'K 

n  im 

R  F  A 
R  I-  A 
OP 
H  F  A 
OP 
PFA 

nn 

T(  I 
FPR 

f  r  r 

RFT 

FKD 


RPU 
EMS 

ru* 

10 

0(  5 

11 

D(  5 

1? 

)=- 

MAT 

VAT 

URN 


TTNF      INPMT(N.R.W.T) 
IflN    PC  IS)  ,W(  IS)  »T(  IS) 

R  .W,T 
,  1^0  )N 
I  =  1  t  M 
•  10DRM  ) 
I  *  1  t  N 
,  1 0  1  )  W  (  I  ) 
I  =  1  .  N 

ni.PG(R(  i  ) ) 

(15) 
(02^.17) 
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c 

C  SUBROUTINE  BESI 

C 

C  PURPOSE 

C  COMPUTE  THE  I  BESSEL  EUNCTION  FOR  A  GIVEN  AROUMFNT  AND  ORDER 

C  ^ 

C  USAGE  r- 

C  CALL  BESK  if,N,BI  ,IER) 

C 

C  DESCRIPTION  OF  PARAMETERS 

C  X   -THE  ARGUMENT  OF  THE  I  BESSEL  FUNCTION  DESIRED 

C  N   -THE  ORDER  OF  THE  I  BESSEL  FUNCTION  DESIRFO 

C  BI  -THE  PESUITANT  I  BESSEL  FUNCTION 

C  IER-RFSULTANT  ERROR  CODE  WHERE 

C  IER=0  NO  ERROR 

C  IFR=1  N  IS  NEGATIVE 

C  IER=2  X  IS  NEGATIVE  

C  IER=3  UNDERFLOW,  BI  .LT.  l.E-69,  BI  SET  TO  0.0 

C  IER=4  OVERFLOW,  X  .GT.  170  WHERE  X  .GT.  N 

C 

C  REMARKS 

C  N  AND  X  MUST  BE  .GE.  ZERO 

C 

C  SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

C NONE  

C 

C  METHOD 

C  COMPUTES  I  BESSEL  FUNCTION  USING  SERIES  OR  ASYMPTOTIC 

C  APPROXIMATION  DEPENDING  ON  RANGE  OF  ARGUMENTS. 

C 

C 

c 


c 


c 


SUBROUTINE  BFSI(X,N,  BWIER) — 

REAL*8  X,XX,TER*,BI,TOL 
C 
C      CHECK  FOR  ERRORS  IN  N  AND  X  A7JD  EXIT  IF  ANY  ARE  PRESENT 

C  

IER  =  0 

BI=1.C 

IF(N)1S0,  15,  1C 

10  IF(X) 160,20,20    

15  IF(X)  160,17,20 
17  RETURN 
C 
C      DEFINE  TOLERANCE 


20      T0L=1.0D-12 


C 

C      IF  ARGUMENT  GT  12  AND  GT  N,  USE  ASYMPTOTIC  FORM 

C 

IFIX-12. ) AC, 40, 30 

30  IF(X-FLOAT(N) 140,40,110 
C 

C      COMPUTE  FIRST  TERM  OF  SERIES  AND  SET  INITIAL  VALUE  OF  THE  SUM 
C 

40  XX=X/2. 

50  TERM=1.0  

IF(N)  70,70,55 

55  DO  6C  1=1, N 
FI  =  I 

I F( DABS (TERM) -1.0-68)56 ,60,60 

56  IER=3 
BI-0.0 
RETURN 

60  TERM=TERM*XX/F  I    ~  

70  BI=TERM 
XX=XX*XX 
C 

C      COMPUTE  TERMS,  STOPPING  WhEN  ABS(TERM)  LE  ABS<SUM  OF  TFRMS) 
C      TIMES  TOLERANCE 


DO  90  K=l,1000  
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c 
c 
c 


IF  (l)ABS(  TCP  M) -OARS  (RT*T0l  )) 100,10^,80 

80rK=K*(N4K)  -       

TFRV=TFRM*{XX/FK) 
9C  HI=MI+TERM 

RETURN  HI  AS  ANSWER 

100  RETURN 

X  f,T  12-AND-*  GT  NrSO-tlSF— ASYMPTDTIC-APPRDX  [MATICN 

I1C  FN=4*N*N 

IF(X-170. 0)115, 111,111 

111  IFR=A 

RETURN 
115  *X=1./(8.*X> 

TERM=1.  

00  130  K=l,30 

IF(l)ABS(TEPM)-DABS(TnL*BI))lA0,lAO,120 
,120  FK=(?*K-1)**? 

TFRM=TFRM*XX*(FK-FN) /FLOAT(K) 
130  8I=BI+TERM  ' 

SIGNIFICANCE  LOST  AFTER  30  TERMS,  TRY  SERIFS  | 

GO  TO  40 
140  PI=3.1415926S3 

BI=BI*DEXP(X)/DSQRT(2.*PI*X) 
GO  TO  100 

150  IFR=1 

GO  TO  100 
160  IER=2 
G  0  -T  0-10  0 1 

END 


\ 
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c 

r. 
C 
C 

c 
c. 
c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

■c- 

c 

c 

c 

c 

c 

c 


SUBROUTINE  PESJ 

PURPOSE 

COMPUTE  THE  J  BESSFL  FUNCTION  FOR  A  GIVEN  ARGUMENT  AND  PROF 

USAGE 

CALL  BESJ( X,N,BJ,D, TER ) 

DESCRIPTION  OF  PAPAMFTFRS 

.  X   -THF  ARGUMFNT  OF  THF  J  BFSSFL  FUNCTION  DESIRED 
N   -THE  ORDER  OF  THE  J  BESSFL  FUNCTION  DESIPFD 
BJ  -THF  RESULTANT  J  BFSSEL  FUNCTION 
D   -REQUIRED  ACCURACY 
IER-RESULTANT  ERROR  CODE  WHERE 
IER=0   MO  ERROR 

-■: IER-1 — N  IS  NEGATIVE — 

IER=2   X  IS  NFGATIVF  OR  7ERO 

IPR=3   REQUIRED  ACCURACY  NOT  OBTAINED 

IER  =  4   RANGE  OF  N  COMPARED  TO  X  NOT  CORRECT  (SFF  PEMAPKS 

REMARKS 

N  MUST  BE  GREATER  THAN  PR  EQUAL  TO  ZERO,  BUT  IT  MUST  BE 

LESS  THAN 

20+10*X--X**— 2A3 FfVR— X-IESS— T-HAN-OR-  EQUAL— TO"  15 

90+X/2  FOR  X  GREATER  THAN  15. 

SUBROUTINES  AND  FUNCTION  SUBPROGRAMS  REQUIRED 

NONE 

METHOD 

RECURRENCE  RELATION  TECHNIQUE  DESCRIBFD  BY  H.  GOLDSTEIN  AND 

RvM.-THALERv/?RECUR-RENCE-TECHNIQ»JES~FOR    THF -CALCULATION  -OF  — 

BESSFL  FUNCTIONS5UM.T.A.C. ,V. 13, PP. 102-l^fl  AND  I. A.  STEGHN 
AND  M.  ABRAM0WITZ»3GENERATI0N  OF  BESSFL  FUNCTIONS  ON  HIGH 
SPEED  C  OMPUT  ER  S3,  M.T.  A.  CV.l  1,1 957,  PP.  2  5  5-2  57 


SUBROUTINE 
REAL* 8 


BESJ<  X,N,BJ,D, IER) 

F  M 1 ,  F  M ,  8  M  K-,-B-d-rAl-F-Anrtrt>R-e\hrD- 


10 

■?,■> 
30 

31 
32 


RJ  =  .0 

IF (N)  10,20,20 

IER=1  - 

RETURN 

[F(X  )  30, 30, 31 

IER=2 

"RETURN 

IF(X-15. )32,32, 


34 


NTEST=20.+10.*X-X**  2/3 

GO  TO  3ft 
34  NTEST=OG.+X/2. 
36  IF(N-NTEST)40,3B,38 

30  IFR=?4 

RETURN 
-40-  lER-C- 


C 

c 

C 


M  =  N+1   , 
8PREV=.0 

COMPUTE  STARTING  VAtUE-OF~M- 

IF( X-5. ) 5C ,60,60 


50  ma  = 

GC- 

60  MA  = 


X+6. 
TO -70 — 


C 

c 

C 


1.4*X+60. /X 
70  MR=N+IF  IX(X)/4+2 
MZERO=MAX0(MA,MB) 

SET  UPPFR  LIMIT  OF  M 

VVAV=NTEST 


6Q 


110     Of)     1QO     M=M7FRfl.MMAXj3 

SET     F(M)  ,F(M-1  ) 


F"l=l. 06-28 

F  M  =  .  C 

ALPHA=. 0 

IF (N»-(M/2)*2) 120 , 110,120 
110    JT=+1 

'  GO    TCTT3C 

120    JT=1 
130    v2=m-2 


00     lftf>    K=1,M2 

VK=M-K  •  . 

RVK=2.*FLCAT(  VK  )*FM1/X-FM 

F  V  =  F  M  1 

FVl  =  8f^K  

T  F  (  ,vk-n-T)  150V1  40  t~T5Tn~ 

140      'J=RMK 
150     JT=-JT 

S=UJT 
160     AL PHA=ALPHA+PMK*S 

HMK=2.*FM1/X-FM 

IF(N)  180,170,180 
17^  *J=BMK 
ISO  ALPHA=ALPHA+BMK 

RJ=RJ/ALPHA 

IF(OABS(BJ-BPREV)-DABS(0*BJ> 1200,200,190 
ISO  BPRFV=BJ 

IER  =  3 
200  RETURN 

ENO 
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APPENDIX  B 

ROOTS  OF  THE  SHIFTED  LEGENDRE  POLYNOMIALS 
AND  CORRESPONDING  WEIGHTS 


* 


The  roots,  x.,  of  the  shifted  Legendre  polynomial  P   (x)  are 

computed  from  the  roots  of  Legendre  polynomial  PN(r)  using  the  formula 

1+r. 
1 
x„  =  —x —  . 
1    2 

The  coefficients  (or  weights)  of  the  shifted  Legendre 

polynomial  P   (x)  in  the  interval  (0,1)  are  then 


A. 

l 


* 


A, 


l 


where  A.'s  are  the  weights  of  the  Legendre  polynomial  PM(r)  in  the 
interval  (-1,1)  and  were  obtained  from  Stroud  and  Secrest  [20]. 
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ROOTS  AND  WEIGHTS 


ROOTS 


0.78867513459481260D  00 
0.21132486540518690D  00 


0.88729833462074120D  00 
0.50000000000000000D  00 
0.11270166537925830D  00 


0.93056815579702600D  00 
0.6o999052179242780D  00 
0.33000947820757180D  00 
0.69431844202973800D-01 


0.95308992296933150D  00 
0.76923465505284130D  00 
0.50000000000000000D  00 
0.23076534494715830D  00 
0 . 46910077030668140D-01 


0.96623475710157570D  00 
0.83060469323313190D  00 
0.61930959304159840D  00 
0.38069040695840140D  00 
0.16939530676686740D  00 
0.33765242898424040D-01 


0.97455395617137920D  00 
0.87076559279969690D  00 
0.70292257568869830D  00 
0.50000000000000000D  00 
0.2970774243H30150D  00 

0.12923440720030290D  00 

0.25446043828620770D-01 


N  =  2 


N  =  3 


N  =  4 


N  =  5 


N  =  6 


N  =  7 


WEIGHTS 


0 . 50000000000000000D 
0 . 50000000000000000D 


00 
00 


0.27777777777777770D  00 
0.44444444444444420D  00 
0.27777777777777770D  00 


0.17392742256872660D  00 
0.32607257743127270D  00 
0.32607257743127270D  00 
0.17392742256872660D  00 


0.11846344252809440D  00 
0.23931433524968290D  00 
0.28444444444444430D  00 
0.23931433524968290D  00 
0.11846344252809440D  00 


0.85662246169584800D-01 
0.18038078652406920D  00 
0.23395696728634530D  00 
0.23395b96728634530D  00 
0.1803807»652406920D  00 
0.85662246189584800D-01 


0.64742483084434500D-01 
0.13985269574463«OOD  00 
0.19091502525255910D  00 
0.20897959183673430D  00 

0.19091502525255910D  00 

0.13985269574463800D  00 
0.64742483084434500D-01 
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ROOTS  AND  WEIGHTS  (continued) 


ROOTS 


0.98014492824876dOOD  00 
0.898333238706dl300D  00 
0.76276620495816400D  00 
0.59171732124782460D  00 
0.40828267875217520D  00 
0.23723379504183530D  00 
0.1O166676129318d50D  00 
0.19855071751231980D-01 


0.98408011975381270D  00 
0.918O1555366331740D  00 
0.dO66857l635O2952OD  00 
0.66212671170190430D  00 
0.50000000000000000D  00 
0.33787328829809550D  00 
0.19331428364970460D  00 
0 . 81984446336681900D-01 
0.15919880246187050D-01 


0.98695326425858540D  00 
.;.93253168334449190D  00 
0.83970478414951200D  00 
0.71669769706462350D  00 
0.57443716949081550D  00 
0.4255628305091O410D  00 
0.28330230293537630D  00 
0.l60295215850487d0D  00 
0. 674683166 55507700D-01 
0.13046735741414190D-01 

0.98911432907302820D  00 
0.94353129988404730D  00 
0.86507600278702430D  00 
0.75954806460340560D  00 
0.63477157797617220D  00 
0.50000000000000000D  00 
0.36522842202382730D  00 
0.24045193539659390D  00 
0.13492399721297520D  00 
0.56468700115952460D-01 
0.10885670926971520D-01 


N  «  8 


N  =  9 


N  =  10 


N  =  11 


WEIGHTS 


0.50614268145188080D-01 
0.11119051722668690D  00 
0.15685332293894350D  00 
0.18134189168918070D  00 
0.1bl34l891689l8070D  00 
0.15685332293894350D  00 
0.11119051722668690D  00 
0.50614268145188080D-01 


0.40637194180787150D-01 
0.903240d0347428300D-01 
0.13030534820146730D  00 
0.15617353852000140D  00 
0.165H967750062960D  00 
0.15617353852000140D  00 
0.13030534820146730D  00 
0.90324080347428300D-01 
O.4O6371941bO787150D-Ol 


0.33335672154344030D-01 
0.74725674575290000D-01 
0.10954318125799080D  00 
0.13463335965499780D  00 
0.14776211235737600D  00 
0.14776211235737600D  00 
0.13463335965499780D  00 
0.10S54318125799080D  00 
0.74725674575290000D-01 
O.33335672154344030D-O1 

0.278342835580d6790L-01 
0.62790184732451900D-01 
0. 9314510546 3866800D-01 
0.11659688229599510D  00 
0.13140227225512.i90D  00 
0.13646254338894990D  00 
0.13140227225512290D  00 
0.1I659688229599510D  00 
0.93145105463866800D-01 
0.62790184732451900D-01 
0 . 2783428355  3086790D-01 
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ROOTS  AND  WEIGHTS  (continued) 


HOOTS 


0.99078031712335930D  00 
0.95205b626lb52V700D  00 
0.88495133709715200D  00 
0.79365897714330640D  00 
0.68391574949908970D  00 
0.56261670423573420D  00 
0.43738329574426560D  00 
0.31608425050090960D  00 
0.206341022d56o911oD  00 
0.11504866290284750D  00 
0.4794137l8l4762ol0D-01 
0.92196828766404250D-02 


0.9920915273592936OD  00 
0.958799199&H48860D  00 
0.90078904536b65450D  00 
0.82117466972016980D  00 
0.72424637-J51822310D  00 
0.61522915797756730D  00 
0.50000000000000000D  00 
0.38477084202243230D  00 
0.27575362448177640D  00 
0.17882-333027982970D  00 
0.99210954633345100D-01 
0.41200800388511100D-01 
0.79084726407060300D-02 


0.99314190434840570D  00 
0.96421744183178630D  00 
0.91360065753488210D  00 
0.84364645240584240D  00 
0.75762431817907670D  00 
0.65955618446394460D  00 
0.55402747435367150D  00 
0.44597252564632830D  00 
0.34044381553605520D  00 
0.24237568182092280D  00 
0.15635354759415710D  00 
0.86399342465117200D-01 
0.35782558168213250D-01 
0.68580956515938870D-02 


N  =  12 


N  =  13 


N  =  14 


WEIGHTS 


0.23587668193255860D-01 
0.53469662997659180D-01 
0.80039164271672700D-01 

00 
00 
00 
00 
00 
00 


0.10158371336153270D 

0.11674626826917710D 

0.12457352290670110D 

0.12457352290670110D 

0.11674626826917710D 

0.1015837133615  3270D 

0 . 80039164271672700D-01 

0.534696b29976591bOD-01 

O.23587668193255860D-01 


0.2024:^002382657900D-01 
0.46060749918864170D-01 
0.69436755109893600D-01 
0.b9072990380972600D-01 
0.10390802376844400D  00 
0.11314159013144830D  00 
0.11627577661543650D  00 
0.11314159013144830D  00 
0.10390802376844400D  00 
0 . 89072990380972600D-01 
0.69436755109893600D-01 
0.46060749918864170D-01 
0.2O242002382657900D-01 


0.17559730165875880D-01 
0 . 40G79043579880060D-01 
0. 60759285 343951540L-01 
0.78601583579096500L-01 
0.92769198738968500D-01 
0.10259923186064770D  00 
0.10763192673157860D  00 
0.10763192673157860D  00 
0.10259923186064770D  00 
0.9276919873&968500D-01 
0.78601583579096500D-01 
0.60759285343951540D-01 
0 . 40079043579880060D-01 
0.17559730165875880D-01 
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ROOTS  AND  WEIGHTS  (continued) 


ROUPS 


WEIGHTS 


N  =  15 


0.99399625901024220D  00 
0.96863669620035250D  00 

0.92410329170521320D  00 

0.8622C88o56d008460lJ  00 
0.7854860863-4269201)  00 
0.69707567353878130D  00 
0.60059704699871720D  00 
0.50000000000000000D  00 
0.39940295300128260D  00 
G.30292432646121800D  00 
0.21451391369573060D  00 
0.13779H3431991470D  00 
0.75896708294786100D-01 
0.31363303799647090D-01 
0.60037409897574110D-02 


0.1537662099805«590D-01 
0.35183023744O54020D-01 
0.5  3579610233585930D-01 
0. 69785 3389b3076900D-01 
0.83134602908496700D-01 
O.93J8O500007780800D-01 
0.99215742663555600D-01 
0.1012d912096278020D  00 
0. 9921574266 3555600D-01 
0.930d0500007780800D-01 
0.83134602908496700D-01 
0. 69785 338963076900D-01 
0.5  3579610233585930D-01 
0.35183023744054020D-01 
0.1537662099005859UD--1 
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APPENDIX  C 

Legendre  Coefficients 

This  table  gives  the  coefficients  of  the  first  fifteen 
shifted  Legendre  polynomials,  P   (x) ,  in  order  of  lowest  to  highest 
power  such  that 

00 

pXT  (x)  =  3~~  a7  xK  , 

N       ^- —  Nk    ' 

K=0 

where  P   (x)  =  P  (l-2x).   The  recurrence  relation  for  these 

coefficients  is 

(N-1)C„  =  (2N-3)C..  ,    -  (N-2)CM  „   -  (4N-6)CM  .    .  . 
Nk         N-1,k        N-2,k         N-1,k-1 

All  of  these  coefficients  are  integers  but  due  to  machine  conversion 
from  hexadecimal  to  decimal,  noninteger  numbers  may  be  printed  in 
this  table.   In  other  words 

-0.13999999999999990D  03  represents  the  integer  -14.   The  coef- 
ficients are  printed  in  order  of  lowest  to  highest  power  of  x.   For 
example 

P   (x)  =  l-6x-6x  . 
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LEGENDRE  COEFFICIENTS 


N  =  1 

0.99999999999999900D  00 
-0.20000000000000000D  01 

N  =  2 

0.99999999999999900D  00 

-0.59999999999999990D  01 

0.59999999999999990D  01 

N  =  3 

0.99999999999999900D  00 
-0.11999999999999990D  02 

0.29999999999999980D  02 
-0.20000000000000000D  02 

N  =  4 

0.99999999999999900D  00 
-0.20000000000000000D  02 

0.69999999999999990D  02 
-0.13999999999999990D  03 

0.70000000000000000D  02 

N  =  5 

0.99999999999999900D  00 

-0.29999999999999980D  02 

0.21000000000000000D  03 

-0.55999999999999980D  03 

0.62999999999999990D  03 

-0.25199999999999980D  03 

N  =  6 

0.99999999999999900D  00 
-0.41999999999999990I>  02 
0.41999999999999990D  03 
•0.16799999999999990D  04 
0.31*999999999999900  04 
•0.27719999999999990D  04 
0.92399999999999990D  03 


N  =  7 

0.99999999999999900D  00 
-0.55999999999999980D  02 

0.75600000000000000D  03 
-0.41999999999999990D  04 

0.11549999999999990D  0^ 
-0.16631999999999990D  05 

0.12011999999999990D  05 
-0.34319999999999990D  04 

N  =  8 

0.99999999999999900D  00 
-0.72000000000000000D  02 

0.12599999999999990D  04 
-0.923^9999999999990D  04 

0.34649999999999980D  05 
-0.72071999999999980D  05 

0.84083999999999980D  05 
-0.51479999999999990D  05 

0.12869999999999980D  05 

N  =  9 

0.999999999999999900D  00 

-0.899999999999999990D  02 

0.198000000000000000D  04 

-0.184799999999999980D  05 

0.900899999999999980D  05 

-0.2522519999999999901)  06 

0.420419999999999980D  06 

-0.411b39999999999970D  06 

0.218789999999999980D  06 

■0.486199999999999970D  05 

N  =  10 

0.999999999999999900D  00 
■0.109999999999999990D  03 
0.296999999999999990D  04 
•0.343199999999999970D  05 
0.210209999999999980D  06 
■0.756755999999999970D  06 
0.168167999999999990D  07 
•0.233375999999999980D  07 
0.196910999999999980D  07 
■0.923779999999999980D  06 
0.184755999999999990D  06 
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LEGSNDRfi  COEFFICIENTS    (continued) 


N  -   11 

0.99999999999999900D 
-0.131999999999999901) 

0.42b'999999999999SOD 
-G.b0059999999999980D 

0.450449999999999BOD 
-0.2oi601599999999bOD 

0.57177119^99999901) 
-0.105j191999999999QD 

G.124710299999^9bOD 
-0.923779999999999bOD 

0.36798759999999970D 
-0.70543199999999980D 

N  =  12 

0.99999999999999900D 
-0.15599999999999990D 

0.60059999999999980D 
-0.1O009999999999980D 

0.90089999999999980D 
-0.49008959999999990D 

0.17153135999999990D 
-0.39907295999999980D 

0.623551499999999&UD 
-0.64664599999999980D 

0.42678635999999990D 
-0.16224935999999990D 

0 . 27041560000000000D 

N  =  13 

0.99999999999999900D 
-0.18199999999999990D 

0.blo9999999999999CD 
-0.l60159999999999bOD 

0.17016999999999990D 
-0.110270159999999bOD 

0.46558511999999980D 
-0.13302431999999J90D 

0.26189l629999999bOD 
-0.35569529999999980D 

0.32720287599999990D 
-0.19469923199999980L 

0.67603899999999980D 
-0. 104005 99999999980D 


00 
03 
04 
05 
06 
07 
07 
08 
Ob 
07 
07 
06 


00 
03 
04 
06 
06 
07 
08 
08 
08 
08 
08 
08 
07 


00 
03 
04 
06 
07 
08 
08 
09 
09 
09 
09 
09 
08 
08 


N  =  14 

0.99999999999999900D 

-0 . 21000000000000000D 

0.10920000000000000D 

.0.24752000000000000D 

0 . 30630600000000000D 

•0.23279255999999980D 

0.1l6  396279999999bOD 

-0.39907295999999980D 

0.9b026930999999990D 

-0.16360143799999990L 

0.19632172559999980D 

-0.16224935999999990D 

0.87bb5069999999990D 

-0.280816199999999601) 

0.40116599999999990D 

N  =  15 

0.99999999S99999900D 
-0.23999999999999980D 

0.14279999999999980D 
-0.37127999999999980D 

0.529073999999999801) 
-0.46558511999999980D 

0.27159131999999990D 
-0.10974506399999990D 

0.31551705899999980D 
-0.65440575199999970D 

0.98160862799999960D 
-0.1j546208399999980D 

0.79096562999999960L 
-C39314267999999970D 

O.11633813999999980D 
-0.15511751999999990D 


00 
03 
05 
06 
07 
08 
09 
09 
09 
10 
10 
10 
09 
09 
08 


00 
03 
05 
06 
07 
08 
09 
10 
10 
10 
10 

11 

10 
10 
10 
09 
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where 


APPENDIX  D 

Negatives  of  Logarithms  of  Roots  of 
Shifted  Legendre  Polynomials 

The  negative  logarithms  computed  in  these  tables  are 

t±   =  -log  x^^   , 
PN*(xi}  =  °   '  i  =  1,  2,  ...  ,  N; 


N  =  3,  A,  ...  ,  15. 

It  can  be  readily  observed  from  these  tables  that  the  times 

or  t  values  interlace  for  different  N,  making  possible  the  use  of 
more  t  values  for  extrapolation  purposes. 
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NEGATIVES  OF  LOGARITHMS  OF  ROOTS  OF  SHIFTED 
LEGENDRE  POLYNOMIALS 


T  =  -LOG(x) 


N  =  2 

0.237400786l51bl920D  00 
0.15543586830764310D  01 

N  =  3 

0.11957401204924280D  00 
0.69314718055994510D  00 
0.21630110809448010D  01 

N  =  4 

0.7195995 915876 3400D-01 
0.40049171327565440D  00 
0.11086339030929460D  01 
0.26674096665219890D  01 

N  =  5 

0 . 48046021994112870D-01 
0.26235921291155620D  00 
0.69314718055994510D  00 
0.14663539074975080D  01 
0.305952276454C52920D  01 

N  =  6 

0.34348454520252990D-01 
0.1d560129738b333bOD  ^0 
0.479U99O102148920D  00 
0.96576881433932520D  00 
0.17755202021841730D  01 
0.33883233221b8l020D  01 

N  =  7 

0.2577539351044751CD-01 
0.13838246257543370D  00 
0.35250852739216200D  00 
0.69314718O5599451OD  00 
0.12137624862357850D  01 
0.20461274134878730D  01 
0.367H949968461940D  01 


N  =  8 

0 . 20054832283440630L-01 
0.10721418964279870D  00 
0.27080371017260680D  00 
0.5247^b25602874090D  00 
0.89579550440199850D  00 
0.143b7091471258830D  01 
0.22860548602906120D  01 
0.39192958006440780D  01 

N  =  9 

0 . 16047962730479190D-01 
0.85540945516813900D-01 
0.21482113346806630D  00 
0.41229833395649940D  00 
0.69314718055994510D  00 
0.10850843404787260D  01 
0.164343800181081301)  01 
0.25012257285301320D  01 
0.41401866208186200D  01 

N  =  10 

0.13132591978187970D-01 
0.69852151264428bOOD-01 
0.17470489635796030D  00 
0.33310114923617220L  00 
0.55436455322540930D  00 
0.85434267909849090D  00 
0.12612407432151420D  01 
0.183073806478711801)  01 
0.26960971741796480D  01 
0.43392173H2598770D  01 

N  =  11 

0.10945 35 3360154720D-01 
0.58125740495859540D-01 
0.1449379H41644620D  00 
0.2750316744149H30D  00 
0.45449006452666760D  00 
0.69314718055994510D  00 
0.10072323073755050D  01 
0.1425^350622303700D  01 
0.20030436426860710D  01 
0.28740687745821190D  01 
0.45203079484393840D  01 
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NEGATIVES  OF  LOGARITHMS  0?  ROOTS  OF  SHIFTED 
LEGENDRF  POLYNOMIALS  (continued) 


N  =  12 

0.92624472047904000D-02 
0.49126661858812220D-01 
0.12222262181691270D  00 
0.23110140981597780D  00 
0.37992054219297830D  00 
0.57515669219094700D  00 
O.82694536135380250D  00 
0.11517464854216140D  01 
0.15782250281801380D  01 
0.21624000844470620D  01 
0.30377764351420280D  01 
0.46864146371634650D  01 

N  =  13 

0.79399104703031280D- 

G .42073611203765560L- 

0.10448418267570260D 

0.19701943975857120D 

0.32262364680743880D 

0.48576046598152780D 

0.693147180559945101) 

0.95510733743158590D 

0.12882474768486800D 

0.17213457580450050D 

0.23105068410166020D 

0.31892975959097460D 

0.48398206080462550D 


02 

01 

00 

0. 

00 

0. 

00 

0. 

00 

0. 

00 

0. 

00 

0. 

01 

0. 

01 

0. 

01 

0. 

01 

0. 

01 

0. 

0. 

0. 

0. 

0. 

N  =  14 

0.68817204656915940D-02 
0.36438447740720890D-01 
0.90361720362063900D-01 
0.17002176740104690D  00 
0.27756763865621490D  00 
0.41618811793798540D  00 
0.59054100076837180D  00 
0.80749793054938170D  00 
0.10775051727753410D  01 
0.14172663523405600D  01 
0.18556355052547450D  01 
0.24487752135585180D  01 
O.333O2947O63O21520D  01 
0.49823254775913270D  01 

N  =  15 

60218359038299730D-02 
3l8656b3915715830D-01 
7893142602d461700D-01 
14825773409602940D  00 
24145253461422690D  00 
36086130375S01570D  00 
50983104019704670D  00 
69314718055994510D  00 
91778446443415710D  00 
11942722523176470D  01 
15393806770033880D  01 
19820162597733620D  01 
25783819644961590D  01 
34621167382948760D  01 
51153725057554800D  01 
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APPENDIX  E 

Legendre  Derivatives 
The  derivatives  of  the  shifted  Legendre  polynomials  tabulated 
in  this  table  are  evaluated  at  the  roots  and  listed  in  order  of 
smallest  to  largest  root.   For  example, 

P5*   (0.5000)  =  -3.75. 
Recalling  that  the  recurrence  relation  for  the  derivative  of 
the  Legendre  polynomial  is 

(x2-l)  PN'  (x)  =  N  xPN(x)-PN_1(x)   , 
and  the  relationships  between  shifted  Legendre  polynomials  and  Legendre 
polynomials  are 

PN*(X)  =  PN(1"2X)   » 
and 

we  can  determine  the  derivative  of  the  shifted  polynomial  from  the 
formula 

-2x(x-l)^  PN*(x)  =  N[(l-2x)PN*(x)-PN*1(x)]   . 


73 


LEGENDRE  DERIVATIVES 


N  =  3 

-0.59999999999999560D  01 

0.29999999999999980D  01 

-0.59999999999999990D  01 

N  =  4 

0.94332756589370900L  01 

-0.37243057339114480D  01 

0.37243057339114550L  01 

-0.94332756589371720D  01 

N  =  5 

-0.13740666961469280D  02 

0.48517780726008250D  01 

-0.37499999999999980D  01 

0.48517780726008380D  01 

-0.13740666961489610D  02 

N  =  6 

0.18915^86873541660D  02 

-0.62770701859863770D  01 

0.42578674578469820D  01 

-0.42578674578469820D  01 

0.62770701859864400D  01 

-0.18915986873542170D  02 

N  =  7 

■0.2495699^249836780D  02 
0.79712213924980360D  01 
■0.50083119828975910D  01 
0.43749999999999980D  01 
-0.50083H9828975930D  01 
0.79712213924980850D  01 
-0.24956992249837180D  02 


M  =  9 

•0.39632641191048770D  02 

0.12128510598590920D  02 

-0.70151138493662570D  01 

0.53499385187671650D  01 

■0.49^18749999999980D  01 

0.53499385187671720D  01 

-0.701511384936626801)  01 

0.12128510598591260D  02 

■0.39632641191052960D  02 

M  =  10 

0.48266539729458070D  02 
-0.14584224552352420D  02 

0.82353816679841920D  01 
-0.60482620784350310D  01 

0.52615698870016580D  01 
-0.52615698870016690D  01 

0.60482620784350360D  01 
-0.82353816679842400D  01 

0.1458422455235282013  02 
-0.48266539729465570D  02 

N  =  11 

-0.5776424^669799700D  02 

0.172b9087349651640D  02 
-0.95906545604470890D  01 

0.6852748003778435023  01 
-0.57293788755183090D  01 

0.54140624999999990D  01 
-0.57293788755183200D  01 

0.685274800377o4640D  01 
-0.95906545604471760D  01 

0.17289087349652200D  02 
-0.57764248669307150D  02 


N  =  8 

O.318627O4825655890D  02 
■0.99233405513432490D  01 
0.59356616715355440D  01 
■0.47776387675098840D  01 
0.47776387675098840D  01 
■0.59356616715355900D  01 
0.99233405513433910D  01 
■0.31862704825656630D  02 
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LUGiSNDRE  DERIVATIVES    (continued) 


N  =  12 

0.b8l25674969396010D  02 

•0.20242269444593o80D  02 

0.11077b6543b970730D  02 

•0.77531502769622o20D  01 

0.b2947105594040460D  01 

■0.571149392612b7630D  01 

0.57114939261267650D  01 

-0.62947105594040800D  01 

0.77531502769623060D  01 

-0.110776654369706401)  02 

0.20242269444595060D  02 

-0. 66125674969409310]}  02 

N  =  13 

-O.79350759O0681939OD  02 

0.2344325H31591050D  02 
-0.1269445H59465570D  02 

0.67436666274612950D  01 
-0.6941793293762b430D  01 

0.bll03956096434350D  01 
-0.56652343749999990D  01 

0.61103^56098434460D  01 
-0.69417932937626740D  01 

0.87436868274bl3890D  01 
-0.12694451159465950D  02 

0.23443251131591910D  02 
-0.79350759066843100D  02 


N  =  14 

0.91439461254679590D  02 
-0.26691692761957000D  02 

0.144397b9l80353500D  02 
-0.98208762794436260D  01 

0.76617260993701340D  01 
-0.b566j90347580b550D   01 

0.bl32U72119930^50D  01 
-0.bl321072119930300D  01 

0.65863903475806610D  01 
-0.76ol726099370l890D  01 

0.9O2J6762794437690D  01 
-0.14439769180354050D  02 

0.2b691b92761959280D  02 
-0.9143946l2547139oOD  02 

N  =  15 

-0.10439175433549440D  03 

0.305873o445766^590D  02 

-0.16312799416455820D  02 

0.10982495249623420D  02 

-0.8449132195H64820D  01 

0.713285o4764050950D  01 

-0.648205145519371&OD  01 

0.62841796874999980D  01 

-0.64820514551937180D  01 

0.71326584764051460D  01 

-0.844913219511b5310I3  01 

0.109824952496236 /OD  02 

-0.±63127994l645o730D  02 

0.30587364457665940D  02 

-0.104391754335544oOD  03 
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APPENDIX  F 


Coefficients  of  the  Polynomials 
PN*(x)/(x-Xi) 
The  coefficients  of  the  polynomials  P   (x)/(x-x.)  are 
tabulated  in  order  of  smallest  to  largest  root,  x  ,  and  lowest  to 
highest  power.   For  example, 

P5*(x)/(x-CL5000)  -  -1.05x°+30.4x1-188.5x2+389„8x3-252x!\ 
These  are  the  coefficients  of  the  polynomials 

rN*U>  £1    K 

D.,X     , 


x-x .    *—z        ik 
i    K=0 

where  x.  is  a  root  of  the  shifted  Legendre  polynomial  P   (x) ,  and 
the  coefficients  were  obtained  by  synthetic  division 

bi,N-l  =  CNN 

bi,N-2  =  CN,N-l+Xibi,N-l 

« 

i,k-l    N,k  i  ik 

bi,0   =  CN1+Xibil   ,   i  =  1,2,  ...,  N. 

The  Legendre  coefficients  C   are  defined  by 

N 
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COEFFICIENTS  OF  THE  POLYNOMIALS 


Pn  (x)/(x-x.) 


N  ='3 

I  =  1 

■0.11270166537925780D  01 
0.122540333075851601)  02 
-0.20000000000000000D  02 

I  =  2 

-0.20000000000000000D  01 
0.20000000000000000D  02 
■0.20000000000000000D  02 

I  =  3 

■0.8d729833462074170D  01 
0.27745966692414800D  02 
■0.20000000000000000D  02 


N  =  ^ 

I  =  1 

-0.10492189413613050D  01 

0.30375707853927490D   02 

•0.1d846520964827400D  03 

0.38982133941172620D  03 

-0.25199999999999980D  03 

I  =  2 

-0.12999934329938890D  01 
0.37309820063988910D  02 

-0.22449609985934440D  03 
0.43615286692668340D  03 

-0.25199999999999980D  03 

I  =  3 


N  =  4 

I  =  1 

-0.10746123148212680D  01 

0.20337454669259780D  02 

-0.74860229094208130D  02 

0.70000000000000000D  02 

I  =  2 

-0.14925584280277500D  01 
0.27623437899478340D  02 
■0.93100663474529990D  02 
0.70000000000000000D  02 

I  =  3 

-0.30302159969205O70D  01 

0.51422110950418290D  02 

•0.116899336525469801)  03 

0.70000000000000000D  02 

I  =  4 

■0.14402613260230380D  02 

0.80616996480843440D  02 

-0.13513977090579180D  03 

0.70000000000000000D  02 


-0.20000000000000000D  01 

0.55999999999999980D  02 

-0.30799999999999980D  03 

0.50400000000000000D  03 

■0.25199999999999980D  03 

I  =  4 

•0.43334063016653720D  01 

0.11122377887464790D  03 

-0.42803749907929230D  03 

0.57184713307331550D  03 

-0.25199999999999980D  03 

I  =  5 

-0.21317381323979400D  02 
0.18509069320743520L  03 

-0.53100119141308810D  03 
0.61817866058827150D  03 

-0.25199999999999980D  03 
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COEFFICIENTS  OF 
* 


THE   POLYNOMIALS 


(x)/(x-xi)    (continued) 


N  =  6 
I  =  1 

-0.10349451752281860D  01 

0.42396585843853390D  02 

-0.39079883163057260D  03 

^'.133425252yl024040L  04 

-0.18791990844381410D  04 

0.92399999999999990D  03 

I  =  2 

-0.12039421497938950D  01 

0.4911b093591293450D  02 

-0.44652276760565730D  03 

0.14850352308906720D  04 

-0.200452126345^-565013  04 

0.9^399999999999990D  03 

I  =  3 

-^.16147012919478840D  01 
0.65210193999593780D  02 

-0.57287955811880220D  03 
0.17876688076475020D  04 

-0.219975793c0295620D  04 
0.92399999999999990D  03 

I  =  4 

-0.26268064067852090D  01 
0.10342575718625110D  03 

-0.53157924924633320D  03 
0.22286370637292510D  04 

-0.24202420639704350D  04 
0.92399999999999990D  03 

I  =  5 

-0.59033512739302930D  01 
0.21309H9724165760D  03 

-0.12214553443508480D  04 
0.27069501770803290D  04 

-0.26154787365474120D  04 
0.92399999999999990D  03 


(N  -  6) 

I  =  6 

-0.29616253702314490D  02 

0.36676017213735010D  03 

-0.15767642490477790D  04 

0.30574561913498280D  04 

-0.27406009155618540D  04 

0.923999999999999990  03 


N  =  7 

1  =  1 

■0.10261104515225810D  01 

0.56409282626533240D  02 

-0.71785734688499980D  03 

0.35730629700534000D  04 

-0.O1852184575646220D  04 

0.66673308224198250D  04 

-0.34319999999999990D  04 

I  =  2 

-0.11484146919320910D  01 
0.62992366443543510D  02 

-0.79586014799722290L  03 
0.39093642194310200D  04 

-0.87746183861063030D  04 
0.90235324855H4300D  04 

-0.34319999999999990D  04 

I  =  3 

-0.14226317870359660D  01 
0.77643498872533810D  02 

-0.96505152144649300D  03 
0.4b021405350141410D  04 

-0.988424572674819301)  04 
0.95995697202363800L  04 

-0.34319999999999990D  04 
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COEtftflCLdJNrS    OJ    THE   POLYNOMIALS 
P      (x)/(x-x.)    (continued) 


(N  =  7) 

I  =  4 

-O.20000000000O00000D  01 

0.10799999999999990L  02 

■0.12959999999999990D  03 

0.58080000000GOOOOOD  04 

-0.11483999999999980D  05 

0.10295999999999990D  05 

-0.34319999999999990D  04 

I  =  5 

-0.336bl^5«5b66?4900D  01 

0.17717224478^888801)  03 

-0.19484070745293940D  04 

0.75791451696147850D  04 

-0.133663971255o6240D  05 

0.10992430279763600D  05 

-0.34319999999999990D  04 

I  =  6 

-0.77378774094586310D  01 

0.37344o388l2587250D  03 

-0.^9601529512276080D  04 

0.95937844698797980D  04 

-0.1513695595o549090D  05 

0.115684675144oo540D  05 

■0.34319999999999990D  04 

I  =  7 

-J.39298839801385420D  02 

o.65633b2l914262670D  03 

-G.39166709579142760D  04 

0.1113450263600b8lOD  05 

-0.16328564345465460D  05 

0.11924669177580150D  05 

•0.34319.999999999990D  04 


M  =  8 

1  =  1' 

■0.10202572815328950D  01 

0.724175993496054«OD  02 

-0.12116395916798290D   04 

0.81909931653316720D  04 

■0.26994994385108750D  05 

0.45990143208138260L  05 

-0.3ti865534773438350D  05 

0.12o69999999999980D  05 

I  =  2 

-0.11131726590019590D   01 

0.78909278079304370D   02 

-0.1314757O994415500D  04 

0.88221628223031720D   04 

-0.28750d421d7334650D   05 

0.482^3928433314890L  05 

-0.39918451217843300D  05 

0.12b69999999999980D  05 

I  =  3 

-0.13H017705ob97280L  01 

0.92674507384851140D  02 
■0.15303843891185150D  04 

0.10107442569908190D  05 
-0.32175727333695730D   05 

0.5^304719856449930L  05 
-0.416b3198942l88390L  05 

0.126b9999999999980D  05 

I  =  4 

-0.16899961587937520L  01 
0.118b236364l6346oOD  03 

-„.19265836709615300D  04 
0.123562655H410070D  05 

-0.37o76325b508b3090I)  05 
0.58128557529128660L  05 

-0.43864598075540460D  05 
0.12669999999999980D  05 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
P   (x)/(x-x- )  (continued) 

Li  -J- 


(N  =  8) 

I  =  5 

-0.2449283430431910UD  OX 
0.170349417o6849940D  03 
■0.2668863H62648720D  04 
0.1609457668843U80D  05 
■0.45447490861672100E  05 
0.652U9690758b5700D  05 
-0.46225401924459470D  05 
O.12869999999999980D  05 

I  =  6 

-o.421525103463<-bi980D  01 

0.2o572973320b59910D  03 
-0.41O679375O1045800D  04 

0.21u37752955b36010D  05 
-0.54&49887815727500D  05 

0.7259552o2C33191yOD  05 
-0.48426o010578ll540D  05 

0.12869999999999980D  05 

I  =  7 

■0.983605641b834bOOOD  01 

0.61144O05&28160780D  03 
-0.63791935^89591110D  04 

0.2bl39054049247060D  05 
-^.6404203171159O04OD  05 

0.78983221126254860D  05 
-0.501715487b2156660D  05 

0.12b69999999999980D  C5 

I  =  8 

-O.503649653H^9uO2OD  02 

0.10896477716111730E  04 
-0.85797840734699980D  04 

0.33251752237512180D  05 
-0.7O4^270OO54007100D  05 

0.83^65934567508040D  05 
-0.51224465226561610D  05 

0.12869999999999980L  05 


=  9 


I  =  1 


■0.10161774228702810D   01 

0.90423351504540900D  02 

■0.19201451290045090D  04 

0.16827750646093990D  05 

.0.744474437429280901)  05 

0.18068097575383680D  06 

■0.24361738382250640D  06 

0.17094402457756930D  06 

-0.48619999999999970D  05 

I  =  2 

-0.108930616262694DOD  01 

0.96850966721181560D  02 
-0.20513258471102790D  04 

^.17895855998660920D  05 
-0.78641525966798870D  05 

0.1891149592o230540D  C6 
-0.25196200632405130D  06 

0.174156O8378088930D  06 
-0.48619999999999970D  05 

I  =  3 

-0.12396401469971320)  01 

0.11003090553602840D  03 

-0.23180887631484420D  04 

0.20034954021466020D  05 

-0.8b643047494983990D  05 

0.20504757819858140D  Ou 

-0.26698430062184840D  06 

0.17956894047104830D  06 

-0.48619999999999970D  05 

I  =  4 

-0.15102o49384090920D  01 

0.133644683861402301)  03 
-0.27885226249108700D  04 

0.23o98b01940943240D  05 
-0.10026992852834280D  06 

0.22953623345146140D  06 
-0.288268877604557101)  06 

0.18659739927705310L  06 
-0.48619999999999970D  05 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
Pn*(x)/(x-xi)  (continued) 


(N  =  9) 

I  =  5 

•0.20000000000000000D  01 

0.17599999999999990D  03 

-0.36079999999999980D  04 

0.29743999999999990D  05 

■0.12069199999999980D  06 

0.263H999999999980D  06 

-0.3H59999999999980D  06 

0.19447999999999970D  06 

-0.48619999999999970D  05 

I  =  6 

-0.29596894298360320D  01 

0.257ol2287l6420720D  03 
-0.50977327077604850D  04 

0.39607355052089010D  05 
-0.14941295064252480D  06 

0.30437164735776020D  05 
-0.3434670826651^310D  06 

0.20236260072^94640D  06 
-0.48619999999999970D  05 

I  =  7 

-0.51729234959793730D  01 

0.438d0397714269750D  03 
-0.7972489H85476270D  04 

0.54354549922924990D  05 
-0.18485674934310280D  06 

0.34863047564048830D  06 
-0.37136171732450760D  06 

0.20939105952895120D  06 
-0.48619999999999970D  05 

I  =  8 

-0.12197435546412570D  02 

0.94899176526824240D  03 

-0.1257565o490961220D  05 

0.72017849395379b70D  05 

-0.22043389218491320D  06 

0.38809931928332460D  06 

-0.39422941985782490D  06 

0.21480391621911020D  06 

-0.48619999999999970D  05 


(N  =  9) 

I  =  9 

•0.62814542856847930D  02 

0.17076420628015200D  04 

■0.17108039318556490D  05 

0.86179083022441700D  05 

-0.24566246209640380D  06 

0.413918811O5224O00D  06 

-0.40836921177951900IJ  06 

0.21801597542243020D  06 

-0.48619999999999970D  05 


(N  =  10) 
I  =  1 

■0.1013219203018738OD  01 

0.11042749919759750D   03 

-0.28973737707332620D  04 

0.31O38008310223160D  05 

•0.18072992729171640D   06 

0.58364067840740420D  06 

-0.11125545264977240D  07 

0.12373488367958950D  07 

-0.74143446270864040D  06 

0.18475599999999990D  06 

I  =  2 

-J.10723496238633440D  01 

0.1168085249-375510D  03 
-0.30596188054902030D  04 

0.33522058019943660D  05 
-0.18947124814716320D  06 

0.6083275903487693'JD  06 
-0.11510090529060490L  07 

0.12683225333985980D  07 
-0.75148917631200480D  06 

0.18475599999999990D  06 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
P   (x)/(x-xi)  (continued) 


(N  =  10) 

I  =  3 

■0.1190894727308265013  01 
0.12958018976025950L  03 
■0.33826409755621850D  04 
0.368431377412866001)  05 
■0.20646168216643710D  06 
0.65534260161554700D  06 
-U.12222597962496650L  07 
0.13236797321288390D  07 
■0.76863950289967240D  06 
0.18475599999999990D  06 

I  =  4 

-0.13952bb4236969410D  01 

0.151534o9oo2067490D  03 

-0.39325717310421160D  04 

J.42399226889406180D  05 

-0.23414442909178530D  06 

0.7291v387497500370D  06 

-0.13289928639733150D  07 

0.14019399534028090D  07 

-0.79136580028112810D  06 

0.18475599999999990D  06 

I  =  5 

-0.17408344256148520D  01 

0.18846128231978200D  03 
-0.48421983559068590D  04 

0.51315971893362j30D  05 
-0.27660619415199490D  06 

0.83585782980166610D  06 
-0.14724363518260360D  07 

0.14994218583338540L  07 
-0.8l764928o31355450D  06 

0.16475599999999990D  06 


(N  =  10) 

I  =  6 

■0.23498292b091?5l60D  01 

0.25295952325131730D  03 

■0.63845812696981860D  04 

0.65643464907114140L  05 

■0.3397'-66772j184200D  06 

0.97999471029732480D  06 

-0.16488406397314120D  07 

0.16094435678254130D  07 

■0.8451547136o644480D  06 

0.18475599999999990L  06 

I  =  7 

-0.35297983448739930D  01 

0.37581834158055780D  03 

-0.91569381241888310L  04 

0.88820534161174900D  05 

-0.42848033560290120D  06 

0.11587468968509620D  07 

-0.18458484019747700D  07 

0.172222955H537820L  07 

-0.87143819971887110D  06 

0.18475599999999990D  06 

I  =  8 

-0.623O4893690947910D  01 

0.64731508099210330D  03 
-0.14490045174987220D  05 

0.12370696236546900D  06 
-0.53963580369c67550D  06 

0.13545020364415510D  07 
-0.20410962474613890D  07 

0.18257797089314560D  07 
-0.89416449710032710D  06 

0.18475599999999990D  06 
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COEFFICIENTS  OP  THE  POLYNOMIALS 
P   (x)/(x-x.)  (continued) 


,1  =  9 

■0.14821771900816640D  02 
0.14107099868099880D  04 
■O.2311144O90272363OD  05 
0.166l3070627bo5280L  06 
■0.65333323711654740D  06 
0.15329086008107720D  07 
-0.22050557441480030D  07 
0.1907625122y025560D  07 
.0.9H314823687994bOD  06 
0.18475599999999990D  06 

I  =  10 

-0.76647524700428040D  02 

0.25563846743443450D  04 

-0.31702590889667050D  05 

0.20061792943536710D  06 

-0.73520846553093450D  06 

0.16515651804509890D  07 

.0.23082263752316030D  07 

0.19570891351267680L  07 

-0.92136953729135910D  06 

0.18475599999999990D  06 


n  =  n 

1  =  i 

■0.10110054727960090D  01 
0.13243059035446710D  03 
.0.42033254270432940D  04 

0.56471403690313600D  05 
■0.39831451706792320D  Oo 
0.16375270636813190D  07 
■0.41250893009936650D  07 
0.64470107363448310D  07 
-0.60903164443090320D  07 
0.31821231006133540D  07 
•0.70543199999999980D  06 


I  =  2 

•0.10598482530719400D  01 
0.13877669109972230D  03 
■0.43996667724858990D  04 
0.58991506942434560D  05 
-0.41488660005838970D  06 
0.16990738941449220D  07 
-0.42591465766413190D  07 
0.66163925077269420D  07 
■0.62050273191705860D  07 
0.32142788280601930L  07 
-0.70543199999999980D  06 

-0.1±559677957135300D  01 

0.1512'9l48747941770D  03 
-0.47b42599947134420D  04 

0.63896975326103250D  05 
-0.44684284782901690D  06 

0.18162255652787900D  07 
-0.45099926736514970D  07 

0.69264750230549010D  07 
-0.64093269944861990D  07 

0.32696237052019400D  07 
-0.70543199999999980D  06 

I  =  4 

-0.131657237583O5340D  01 

0.17205419079357020D  03 
-0.54215736977179900D  04 

0.71935442730423050D  05 
-0.49834180996460860D  06 

0.20007610589184740D  07 
-0.48936349314803500D  07 

0.73837131971996860D  07 
-0.66978207698503270D  07 

0.33440664896906870D  07 
-0.70543199999999980D  06 
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COEFFICIENTS  OP  THE  POLYNOMIALS 
P   (x)/(x-x.)  (continued) 


(N  =  ID 


■0.15753698412222740D  01 

0.20546702890291260D  03 

-0.64346500580881540D  04 

0.84479759022740410D  05 

-0.5765384804153b6lOD  06 

0.22708602111336490D  07 

-0.54300663552956530D  07 

0.79900452708907000D  07 

-0.7059208201154o8lOD  07 

0.34320878162051090D  07 

•0.70543199999999980D  06 

T   =   6 

■0.20000000000000000D  01 

0.25999999999999980D  03 
-0.805999999999999SOD  04 

0.10399999999999990D  06 
-0.69289999999999980D  06 

0.26502319999999990D  07 
■0.61349599999999980D  07 

0.87339199999999990D  C7 
-0.74742199999999980D  07 

0.35271599999999980D  07 
-0.70543199999999980D  06 

I  =  7 

-0.273801254146d0460D  01 

0.35392094279590580D  03 
■0.10777033822814850D  05 

0.13493737947363220D  06 
-0.86387751199106730D  06 

0.31600456547536580D  07 
-0.70029225301624410D  07 

0.95802989549627140D  07 
-0.79148578553086700D  07 

0.36222321837948840D  07 
-0.70543199999999980D  06 


(N  =  11) 

•I  «  8 

■0.41588353129710050L  01 

0.53167035015198220D  03 

-0.15630274065580460D  05 

0.18477591316176540D  06 

-0. 11048947740846 370D  07 

0.37975207993617690D  07 

■0.79857589728738290D  07 

0.10464299332738090D  08 

•0.83456623626341160L  07 

0.37102535103093070D  07 

-0.70543199999999980D  06 

1-9 

-0.741158000545663901)  01 

0.92339704254301420D  03 

-0.24951847165799960L  05 

0.26020o88357447880D  06 

-0.14100020778751850D  07 

0.45063438282596450D  07 

-0.89781521209175620D  07 

0.11^93527545564720D  08 

-0.872715364766870601)  07 

0.37846962947980530L  07 

-0.70543199999999980D  06 

I  =  10 

-0.17708925439165380D  02 

0.^u239721177599240D  04 

-0.40128918809659610D  05 

0.35295803072169190D  06 

-0.17264780148669660D  07 

0.51628244414053760D  07 

-0.98264b24022726930D  07 

0.11961628235470690D  08 

-0.90209578666288160D  07 

0.38400411719398010D  07 

-0.70543199999999980D  06 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
P   (x.)/(x-x.)  (continued) 


(taJ.1) 

i  -  11 

-0.91863882962169210D 

0.36870595581192030D 

■0.55388450186096020D 

0.42914670535641570D 

■0.195700336584b80dOD 

0.56048574b30623620D 

-0.10367254135710960D 

0.12370929196046650D 

-0. 91956485 387888340D 

0.38721968993866440D 

-0.70543199999999980D 

N  *  12 


02 
04 
05 
06 
07 
07 
08 
08 
07 
07 
06 


I  -  1 


-0.10093054785081420D 

0.15643295677440690D 

-0.59040000514031970D 

0.95072538604810130D 

-0.81332606983436670D 

0.41256067157586650D 

-0. 1314875665 J783650D 

0.27007540306116740D 

-0.35676536042330500D 

0.29257811703237760D 

-0.1354571146C768950D 

0 .  27041560000000000D 

1  =  2 

-0.10503554820636720D. 

0.16275220866733050D 

-0.61374873546083530L 

0.98694040328690540D 

•0.84260142802385u90D 

0.42626519542308540D 

■0.13539590592588070D 

0.27695463941829490D 

■0.36404991281091360D 

0.29682634957866470D 

-O.1365042O94824111OL 

0 . 270415  60000000000D 


01 

03 

04 
05 

06 
07 
08 
08 
08 
08 
08 
07 


01 
03 
04 
05 
06 
07 
08 
08 
08 
08 
08 
07 


(N  *  12)  -. 

1-3 

-0.11300056386074940D  01 

0.17500396673720190D  03 

-0.65890583909278350D  04 

0.10566789120383710L  06 

-0.89861676621079580D  06 

0.45225980977865030E  07 

-0.14272579036542980D  08 

0.28967374689262430D  08 

■0.37728374331019390D  08 

0.30438086864005020D  08 

-0.13831889532080700D  08 

0.27041560000000000D  07 

I   St  > 

-0.12599870085909810L  01 

0.19497040599021400D  03 
-0.73218217916793060L  04 

0.11689929917036390D  06 
-0.98783069732489310L  06 

0.49304H4429093590D  07 
-0.15400474144556430D  08 

0.30878277145749980D  08 
-0.39660450849491590D  08 

0.31504903076265990D  08 
-0.14078758315004040D  08 

0.27041560000000000L  07 

I  -  5 

-0.14621b8404l825510D  01 

0.22596033460115980D  03 
-0.84513913733266510D  04 

0.13400569981580050D  06 
-0.11213286150317220D  07 

0.55263640115562260D  07 
-0.17000298643450430D  03 

0.33493887768086840D  08 
-0.42200025738625690D  08 

0.32846990697067070D  08 
-0.14375521122497510D  08 

0.270415&OOOOOOOOOOD  07 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
?n  (x)/(x-xi)  (continued) 


(N,  m   12) 

1  =  6 

-0.17774090110084440D  01 

0.2741l6622955oOOOOD  03 
-0.10187901165549110D  05 

0.1598l05746d493100D  06 
-0.13172190226655770D  07 

0.63696597527710040D  07 
-0.19166647853966690D  08 

0.36864bl4913043360D  08 
-0.45307106764057200D  08 

0.34406182911952610D  08 
-0.14703532663486610D  08 

0.270415bOOOOOOOOOOD  07 

1=7 

-0.22863241685604960D  01 

0.35143929209705040D  03 
-0.12928l5880926blOOD  05 

0.19930308733532540D  06 
-0.160407797603712201)  07 

O.75375947230744230D  07 
-0.219d4244415561060D  08 

O.4O9779O60124837OOD  08 
-0.48875309586617870D  08 

O.36099436277086250D  08 
-0.15042183336513350D  08 

C.27041560000000000D  07 

1  =  8 

-0.31b37134669490470D  01 

0.48353O2179429881OD  03 

-0.17471512020309010D  05 

0.26141286017492700D  06 

-0.^0231540762G48540D  07 

0.91043508787125080D  07 

-0.25464049880790580D  08 

0.45694292253791910D  08 

-0.5271018002^588890D  08 

0.37820359472G91660D  08 

■0.15370194877502440D  08 

0.27041560000000000D  07 


(N  =  12) 


•0.48463460447902640D  01 

0.73254291300179080D  03 

■0.25556997895957640D  05 

0.3ol26ll8341389720D  06 

■0.2bl52764443786540D  07 

0.11076903293286250D  08 

-0.294475263454216901)  08 

0.50b9l663294906840D  08 

-0.56525292661719900D  08 

0.39445899926225400D  08 

-0.15666957684995930D  06 

0.27041560000000000D  07 

I  *  10 

-0.86-19741157220160D  01 

0.12803975480242770D  04 
-0.41074814193766260D  05 

0.5130453&720347610D  06 
-0.3371222255o908350D  07 

0.13^95884596250260D  08 
-0.33527129359225730D  Oo 

0.55456243295603750D  08 
-0.59965118501394010D  06 

0.40847771543197920D  00 
-0.159138264b7919280D  08 

0.27041560000000000D  07 

I  =  11 

-0.20858610712861270D  02 

0.26188844668536790D  04 
-0.6647943920880672UD  05 

0.70128491360441830L  06 
-0.41637333025608900D  07 

0.15376337170479350D  08 
-0.37061910459839990D  08 

0.59351358387366690D  08 
-0.62655520668857870D  08 

0.41907005475455270D  08 
-0.16095295051758860D  06 

0.27041560000000000D  07 
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(N 


COEFFICIENTS  OP  THE  POLYNOMIALS 
?n  (x)/(x-x. )  (continued) 

12)  (N  - 


IV 


1-12 

•0.10846360047086500D  03 
0.51559690463514850D  04 
•0.92197417744400200D  05 
0.d571425244595l830L  06 
■0.47460933446364170D  07 
0.16790453363184500L  08 
■0.393378646172721801)  08 
0.61762577991757680D  08 
■0.64272493552206250D  08 
0.42529277095548140D  08 
•0.16200004539231030D  08 
0.27041560000000000D  07 


-N  »  13 

1  =  1 

-0.10079715286689510D  01 

0.l824348091o835120D  03 

■O.807139761o066i8701j  04 

0.15330097897999020D  06 

■0.15607421072745860D  07 

0.95417344384769900D  07 

-0.373H857364665380D  08 

0.964754359bl292820D  08 

-0.16673481171541170D  09 

0.19042647081912740D  09 

-0.13786672036695810D  09 

0.57285552860546860D  08 

-0.10400599999999980D  08 


-0.10429712616983220D  01 

0.18873297851273380D  03 

■0.834509146933938501)  04 

0.15833858496354280D  06 

-0.16096815846966000D  07 

0.98220090495688290D  07 

-0.38315116413652590D  08 

0.98778976478833180D  08 

■0.17012180818179740D  09 

0.19350609793310380D  09 

-0.13944189578075440D  09 

0.57631813044520720D  08 

-0.10400599999999980D  08 


-0.11101378477288740D  01 

0.20U81267981965990D  03 

-0.88690991095794720D  04 

0.16795375306639010D  06 

-0.17026697369630140D  07 

0.10351309566871580D  08 

-0.40194985295797750D  08 

0.10305335658962990D  09 

-0.17633237685047210D  09 

0.19907316154864720D  09 

■0.14224164371270630D  09 

0.58235153454759510D  08 

-0.10400599999999980D  08 

I .«  4 

-0.12177677100408040L  01 

0.22015076567274520D  03 
-0.97054250797130950D  04 

C.18321872369925510D  06 
-0.18491574719641470D  07 

0.11176499795303420D  08 
-0.43087072104590820D  08 

0.10952267673582370D  09 
-0.18554999183803190D  09 

0.20714875218926850D  09 
-0.14619803585958360D  09 

0.59063190730108370D  08 
-0.10400599999999980D  08 


87 


COEFFICIENTS  OP  THE  POLYNOMIALS 
Pn  (x)/(x-x.)  (continued) 


(y  ■  13) 

1  =  5 

-0.13807456076844240D  01 

0.24938924169703470D  03 

-0.10963963406266510D  05 

0.20600177182381960D  06 

-0.20651787550969150D  07 

0.12374017389442320D  08 

-0.47200090695790450D  08 

0.11850142742212430D  09 

-0.19798539202253470D  09 

0.21770203249501490D  09 

-0.15119280842328430D  09 

0.60071303146785120L  08 

-0.104vJ0599999999980D  08 

1  =  6 

-0.16254106083538880D  01 

0.29318277109067540D  03 
-0.12835570496802190D  05 

0.23946269059726430D  06 
-0.23767360347638980D  07 

0.14060256821493990D  08 
-0.52823008723021120D  08 

0.13036006216061530D  09 
-0.21379280571123440D  09 

0.23058480315709960D  09 
-0.15704404056613840D  09 

0.61205147619538480D  08 
-0.10400599999999980D  08 

1=7 

-0.20000000000000000D  01 

0.35999999999999990L  03 

-0.15659999999999980L  05 

0-28899999999999990D  06 

-0.28253999999999980D  07 

0.16403231999999980D  08 

-0.60310559999999990D  08 

0.14542751999999990D  09 

•0.23292d21999999980D  09 

0.24545415999999970D  09 

•0.16349743199999980D  09 

0.62403599999999980D  08 

•0.10400599999999980D  08 


(N  =  13) 

1  =  8 

-0.25989495325629240D  01 

0.46625427624523020D  03 
-0.20073625338025120D  05 

0.36407741793960490D  06 
-0.34764135843287480D  07 

0.19623634618423160D  08 
-0.70002386979225560D  08 

0.16379082336259«OOD  09 
-0.25495904554972060D  09 

0.2617044834296682Q]D  09 
-0.17022701675121470D  09 

0.63602052380461440D  08 
-0.10400599999999980D  08 

1  =  9 

-0.36264255887223260D  01 

0.646858494594567701)  03 
-0.27354641374456070D  05 

0.48160875083736360D  06 
-0.44245701265234570D  07 

0.23943278663642210D  08 
-0.82012460865595080D  08 

0.18499071129263050D  09 
-0.27887545939564370D  09 

0.27843637866464460D  09 
-0.17684807380864750D  09 

0.64735896S53214790D  08 
-0.10400599999999980D  08 

I  ■  10 

-0.55920489476264360D  01 

0.986481897035092301)  03 

-0.40282425626881880D  05 

0.67036l262498bl210D  06 

-0.576729b7016894390D  07 

0.29412607766919960D  08 

-0.95880735722672390D  08 

0.20770874136904660D  09 

-0.30299336534800010L  09 

0.29448811625060740D  09 

-0.18294253782839130D  09 

0.65744009269b91540D  08 

•0.10400599999999980D  08 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
?n  (x)/(x-xi)  (continued) 

13)  N  =  14 


I  «  11 

•0. 100795 32080864670D  02 

0.17328778717481680D  04 

-0.65084769641774720D  05 

0.95831383449131930D  06 

-0.74929«47037156350D  07 

0.356213ol825677550D  08 

■0.11024115446467420D  09 

0.2^964364791696370D  09 

-0.32504456994910b90D  09 

0.30b54l8355b664090D  09 

-0.18809455571035100D  09 

0.66572046545240380D  08 

-0.10400599999999980D  08 

1  =  1? 

-0.24^71373142518990D  02 

0.38282903577150800D  04 
-0.10586468226722040D  06 

0.13178219165339170D  07 
-0.9317248203b267870L  07 

0.41498412167982510D  08 
-0.12281557116129330L  09 

0.24778O3523825O51OD  09 
-0.34250008457189330D  09 

0.31929514242560010D  09 
-O.191931552291O2610J  09 

0.67175386955479200D  08 
-0.10400599999999980D  08 

1-13 

-0.12644666618088360D  03 

0.70245338566609130D  04 
-0.14736930837191630D  06 

0.16173403145188670D  07 
-0.10657000989156780D  08 

0.45522697896197000D  08 
-0.1309752402C902070D  09 

0.25909930830793510D  09 
-0.35307976886615020D  09 

0.32566732552037570L  09 
-0.19416523890094200D  09 

0.67521647139453030D  08 
-0.10400599999999930D  08 


•0.10069054782458360D  01 

0.21043b28670453990D  03 

-0.10783518111968030D  05 

0.^3637125475372330D  06 

-0.28441945032009650D  07 

0.20576174872216570D  08 

■0.96481786447879590D  08 

0.30468070295618890D  09 

■0.66011574395699H0D  09 

0.982b3765909997460D  09 

-0.98735094411646760D  09 

0.63952860422322550L  09 

-0.^4097472348001660D  09 

0.40116599999999990D  08 

1  =  2 

-0.1037110468422d440D  01 

0.21671759964702230D  03 
-O.IHOO486I9325809OD  05 

0.24519315203176540D  06 
-0.29224391985846560D  07 

0.21112267750251590D  08 
-0.98820046304836740D  08 

0.3H39543910837490L  09 
-0.67295388233064620D  09 

0.998d0012110107130D  09 
-0.10002071037698320D   10 

0.64537983781746240D  08 
-0.24213507457301090D  09 

0.40116599999999990D  08 


89 


COEFFICIENTS  OF  THE  POLYNOMIALS 
Pn  (x)/(x-xi)  (continued) 


(N  »  1*0 


(N  »  1*0 


•  1  =  3 

■0.10945701042210040D  01 
0.22866164573420590D  03 
■U.11702419723639010D  05 
0.25811888195511410D  06 
•0.30702047934305350D  07 
0.22120224016802250D  08 
-0.103191755835178001)  09 
0.32386273H8846210D  09 
-0.69659163832994400D  09 
0.10282640822576100D  10 
-0.10233718266634350D  10 
0.65578080356590450D  09 
-0.244l6564786l93blOD  09 
0.40116599999999990D  08 


•0.131991539o2437980D  01 

0.275440055S7036710L  03 

■0.14049918526392390D  05 

0.30816075444191710D  06 

-0.36362339215554540L  07 

0.25927127214786090D  08 

-0.11941162739159840D  09 

0.36912929785643290D  09 

-0.78025480169954270D  09 

0.11295302404722740D  10 

■0.11003963250961380D   10 

0.689124230012447201)  09 

■0.25042288827733710D  09 

0.40116599999999990D  08 


I  =  4 

•0.11853306384700660D  01 

0.24751442830b26710D  03 

-0.12650424287634700D  05 

0.27839810745672320D  06 

-0.33007450983788320D  07 

0.23681141365139410D  08 

-0.10989809578462980D  09 

0.34276783051601440D  09 

-0.73194343166031770D  09 

0.10716230071987210D  10 

-0.1O568339927924O10D  10 

0.67049367136484310D  09 

-0.24697197272741550D  09 

0.40116599999999990D  08 


-0.15161710607407580L  01 

0.3io09714612803210D  03 

•0.16077330637859610D  05 

0.350906677571745701)  06 

•0.41120883805108430L  07 

0.29060704866359880D  08 

■0.13241567161503020D  09 

0.40429806385895040D  09 

-0.84294751415744180L  09 

0.1202425031443H00D  10 

■0.115349H543210250D   10 

0.7110879356d535140D  09 

-0.25435704837033370D  09 

0.40116599999999990D  08 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
Pn  (x)/(x~xi)  (continued) 


(N  =  1*0 

1  =  7 

-0.1b049646403537720D  01 

0.37578467674825420D  06 
-0.19031935799849300D  05 

0.41241287621467630D  06 
-0.47843243277376570D  07 

0.33382697661047380D  08 
-0.149»3658064213530D  09 

0.4498628513b408280D  09 
-0.92126560909521140D  09 

0.12900962569385440D  10 
-0.12149596007792130L  10 

0.73558445760512930D  09 
-0.25859050142234330D  09 

0.40116599999999990D  08 

1  =  8 

-0.22422905954241500D  01 

0.46585315789012260D  03 
-C.23441235145503930D  05 

0.50244970703015500D  06 
-0.57416323780458560D  07 

0.39324448510673690D  08 
-0.17281744290778780D  09 

0.50733061810099180D  09 
-0.10156201690733830D  10 

0.13911041045129520D  10 
-0.12828439390025380D  10 

0.76159324053700740D  09 
-0.26292529857765640D  09 

0.4011b599999999990D  08 


(N  «  14) 


■0.29373422408393690D  01 

O.608213b9H324455OD  03 

-0.30289244915877940D  05 

0.63808107291382380D  06 

-0.71229930356286800D  07 

0.47456473659044190D  08 

•0.20249980523924240D  09 

0.57740263089002750D  09 

-0.11^46104691522120D  10 

0.150216O3094537710L  10 

-0.13542^85457758150D  10 

0.78789815524134710D  09 

-0.26715875162966600D  09 

0.40116599999999990D  08 

I '«  10 

-0.4125b264545699580D  01 

0.84940111152544690D  03 

-0.41549543307381450D  05 

0.84979835908124590D  06 

-0.91315334289765970L  07 

0.58371048055375070L  08 

-O.2394O2O369894O50OD  09 

0.65877451818192600D  09 

-0.124391518799657601)  10 

0.16177332191812960D  10 

-0.14254071787365220D  10 

0.81314437068440010D  09 

-0.2710929H72266260D  09 

0.40116599999999990D  08 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
P   (x)/(x-xi)  (continued) 


(N  «  14) 

I  »  11 

■0.63957614994190500D  01 

0.13022041497201630D  04 

■0. 61513128408473^901)  05 

0.11896555879521140D  07 

■0.11981847811414090D  08 

0.72255528335757830D  08 

•0.28231372004949410D  09 

0.74676425157665550D  09 

■0.13655274325948350D  10 

0.17299700043088120D  10 

-0.14917937922113510D  10 

0.83592479863585380L  09 

-0.27454382727258420D  09 

0.4011o599999999990D  08 


(N  «  14) 

I  *  13 

■0.27946576521975180D  02 

0.50877699303161640D  04 

•0.16299086393618340D  06 

0.23622999693438930D  07 

-0.1958384382O272490D  08 

0.10327411926099290D  09 

-0.36671946922631620D  09 

0.904169^0565976880D  09 

-0.15678030641774040D  10 

0.19062727573007680D  10 

-0.15914037904035320D  10 

0.86885374294133100D  09 

-0.27938072542698890D  09 

0.40116599999999990D  08 


I  =  12 

■0.11574162157586150D  02 

0.22966128234427870D  04 

■0.99808481529113370D  05 

0.17096370673251690D  07 

-0.15664736490572830D  08 

0.88131683554205400D  08 

■0.32713902258233050D  09 

0.832p7505630568550D  09 

•0.14779539988498060D  10 

U.18294164821220520D  10 

-0.15486318536737270D  10 

0.85488782922^b7090D  09 

■0.27735015213806370D  09 

0.40116599999999990D  08 


I  =  14 

■0.14501307272487430D  03 
0.93592930947540420D  04 
■0.22757146947683940D  06 
0.29087565319279450D  07 
■0.22499462811690370D  08 
0.11370404087734690D  09 
■0.392563658984134801)  09 
0.949141182414434701)  09 
-0.16226264769256810L  10 
0.19521312846090240D  10 
■0.16164795526579080D  10 
0.87592672246302460D  09 
-0.28054107651998310D  09 
0.40116599999999990D  08 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
Pn  (x)/(x-xi)  (continued) 


N  =  15 

1  =  1 

-0.10060401524772970D  01 
0.2404374b422705070D  03 
-0.14124361524008780D  05 
0.35931285982064320D  06 
-0.496l212977893b520D  07 
0.41848546858240950D  08 
■0.23113042032021540L  09 
0.87155279693145960D  09 
•0.22974108527756610D  10 
0.42722964284119880D  10 
■0.55772743622880000D  10 
0.49989464172226740D  10 
■0.292829058096b2000D  10 
0.10091951654130530D  10 
■0.15511751999999990D  09 


■0.10323781804559840D  01 
0.24670510910533980D  03 
-0.14487676283525810D  05 
0.36834483466923620D  06 
-0.50817764644264680D  07 
0.42819702885791220D  08 
■0.236178969898214201)  09 
0.88915862209255000D  09 
•0.23393827394690690D  10 
0.43408171474656100D  10 
■0.56525518329132970D  10 
0.50521073445626440D  10 
■0.29500729908815030D  10 
0.10131288779044060D  10 
■0.155H751999999990D  09 


(H.s.15) 

I  -  3  . 

-0.10821295604528690D  01 

0.25854022227177740D  03 

■0.15173043861639020D  05 

0.38535406088777140D  06 

-0.53082658433782880D  07 

0.44638133558104910D  08 

■0.24559287741861270L  09 

0.9^182093736457290D  09 

■0.24167749132397410D  10 

0.44662567959738930D  10 

■0.57892115871097700D  10 

0.51476895013676600D  10 

■0.29888074454704940D  10 

0.10200367891668490D  10 

-0.15511751999999990D  09 

1  =  4 

■0.11598116781339040D  01 
0.27700966416469860D  03 
•0.1o240833159132670D  05 
0.41177860837793960D  06 
-0.56586769004905530D  07 
0.47436110584700210D  08 
-0.25997785262680260D  09 
0.97131080496675650D  09 
-0.25328663064844250D  10 
0.46522268248212260D  10 
■0.59891050309551950D  10 
0.52853821741328200D  10 
-0.30436640474547060D  10 
0.10296376990336900D  10 
■0.15511751999999990D  09 
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COEFFICIENTS  OF   THE  POLYNOMIALS 
P   (x')/(x-x.)  (continued) 


(N  =  15) 

1  =  5 

•0.12730970773184230D  01 

0.30392251000383360D  03 

■0.17792902679660220D  05 

0.45002337213037700D  06 

■0.61627019399481760L  07 

0.51427785627769750D  08 

■0.2802895407200716GD  09 

0.10403253648001450D  10 

-0.26924031654720320D  10 

0.49035297017800650D  10 

•0.62541611670470500D  10 

0.54641925653287200D  10 

■0.31133126064361470D  10 

0.10415367462979740D  10 

■0.155H751999999990D  09 


1  =  6 

-0.14345645273585880D 

0.34223749949777700D 
-0.199946l89912b6290D 

0.50394152936850410D 
-0.68669710511210130D 

0.56940074737338770D 
-0.30793105169337000D 

0.11326167563681-360D 
-0. 2901483885 3350630D 

0.52255067461657440D 
-0.65854823344066340D 

0.56819159222852080D 
-0.31958357783530470D 

0. 105525 ^7502683340D 
-0.15511751999999990D 


01 

-0. 

03 

0. 

05 

-0. 

06 

0. 

07 

-0. 

08 

0. 

09  • 

-0. 

10 

0. 

10 

-0. 

10 

0. 

10 

-0. 

10 

0. 

10 

-0. 

10 

0. 

09 

-0. 

(N  -  15) 


■0.16b50098603607240D  01 
0.39683010652589540D  03 
-0.23115614641881100D  05 
0.57969713154261050D  06 
■0.78439327865484000D  07 
0.64460155784839010D  08 
-0.34487542895894950D  09 
0.12530451403345920D  10 
-0.3l67057ol23719550D  10 
0.56227381145203220D  10 
■0.69819660060510230D  10 
0.59344986988532170D  10 
-0.32886568640604680D  10 
0.10702182755502340D  10 
-0.15511751999999990D  09 

I  »  8 

20000000000000000D  01 
47600000000000000D  03 
27607999999999990D  05 
68734400000000000D  06 
92067919999999990D  07 
74703439999999970D  08 
39377575999999980D  09 
14073497599999980D  10 
34956416599999970D  10 
60968317199999950D  10 
74385091199999960D  10 
62153985599999960D  10 
33885154799999980D  10 
10858226399999980L  10 
15511751999999990D  09 
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COEFFICIENTS   OF    THE   POLYNOMIALS 
P      (x)/(x-x.)    (continued) 


(N  ■  15)- 

1  =  9 

■0.25037371216725420D  01 

0.59462820966565230D  03 

■0.34264573377579380D  05 

0.84379803426575600L  06 

-0.11133973678256590D  08 

0.88b93731620030380D  08 

-0.45792748152110470D  09 

0.16011978721570450D  10 

-0.38907391799828800D  10 

0.66432116239475960D  10 

•O.7944044O59289193OD  10 

0.65151354569540330D  10 

■0.34915136019074060D  10 

0.11014270044497620D  10 

■0.15511751999999990D  09 

I  -  10 

-0.33011544884820390D  01 
0.78137945630398640D  03 
-0.44561031797570630D  05 
0.10785497883883430D  07 
-0.13905090623848970D  08 
0.1077939b854364180U  09 
■0.54072029595592120D  09 
0.18378528741743480D  10 
-0.43486692904945690L  10 
0.72473157080254900D  10 
•0.84799085038268430D  10 
0.68211751770212770D  10 
■0.35932443448646990D  10 
0.11163925^97316650D  10 
■0.15511751999999990D  09 


(N  «_15) 

I  *  11 

-0.4ool7u22773344^50D  01 

0.10970770784429030D  04 

-0.61454861805635260D  05 

0.14443125523029000D  07 

-0.17930899592615290D  08 

0.13345340595477830D  09 

■0.6439578284938^780D  09 

0.21140484721628930D  10 

■0.485340134772725101)  10 

0.78813357285103280D  10 

-0.90192310519957370D  10 

0.71183137806629320D  10 

-0.36890032245624610D  10 

0.11301065337020220D  10 

■0.155H751999999990D  09 

1  =  1? 

■0.72573609683551000D  01 

0.1689097344182373OD  04 

-0.91376725490805910D  05 

0.20313590993403970L  07 

-0.23654503729479240D  08 

0.16622^65564157130D  09 

■0.76469843200354300D  09 

0.24149028864505400D  10 

-0.53723899378805530D  10 

0.85032145783714720D  10 

■0.95279838438689320L  10 

0.73896231506956810D  10 

-0.37740682800167090D  10 

0.11420075809663060D   10 

-0.15511751999999990D  09 
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COEFFICIENTS  OF  THE  POLYNOMIALS 
•Pn  (x)/(x-xi)  (continued) 


(N  «  15) 

I  »  13 

-0.13175801987562450D  02 

0.29885907189998530D  04 
-0.14877337284699690D  06 

0.29317032602886430D  07 
-0.31082200963825030D  08 

0.20391280971058980D  09 
-0.891718b5038657040D  09 

0.27106839576540920D  10 
-0.58564678539087620L  10 

0.9059545289H32640D  10 
-0.99680342913989280D  10 

0.76179075692639350D  10 
-0.38440235063014350D  10 

0.11516084908331470D  10 
-0.15511751999999990D  09 

I  m    14 

-0.3l8d4396056873630D  02 

0.66356403417380920D  04 
-0.24373579094520960D  06 

0.40666700762636210D  07 
-0.390287302497177bOD  08 

0.24008254354781750D  09 
-0.10046383076688750D  10 

0.2959265163^077490D  10 
-0.62463262239H8420D  10 

0.94929825630010400D  10 
-0.10301966880242760L  11 

0.77874933494715830D  10 
-0.38950918981242010D  10 

0.115851640209559O0D  10 
-0.15511751999999990D  09 


(N  «  15) 

.1  *  15 

•0.16656281503583290D  08 

0.1223190425593b6lOD  05 

■0.34113659259376500D  06 

0.50207707923543650D  07 

-0.4496683H98449170D  08 

0.26511483494478620D  09 

-0.10787415823338880L  10 

0.31162333H5351080D  10 

-0.64855027109461630D  10 

0.97530538298919200D  10 

-0.10498862328606550D  11 

0.78860949321775700D  10 

-0.39244477505985050D  10 

0.11624501145869460D  10 

-0.15511751999999990D  09 
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APPENDIX  G 


ELEMENTS  OF  THE  EXPLICIT  INVERSE 
MATRIX  WITH  DIVISION  BY  WEIGHTS 


From  equations  (37j  and  (45)  it  is  evident  that  the  explicit 
Inverse  matrix  with  division  by  weights  Is  the  necessary  matrix  for 


use  in  obtaining  the  temperature  function; 


a(x, )  =  u(-log  xj 
i        6  i 


N-l 


<-0 


F(k+1) 


q 


<1 


A 
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ELEMENTS  OF  THE  EXPLICIT  INVERSE 
MATRIX  WITH  DIVISION  BY  WEIGHTS 


3 


;N 


1  =  1 


1  =  1 


0.67620999227555400D  00 

•0.73524199845511530D  01 

0.12000000000000080D  02 

1  =  2 

-0.14999999999999990D  01 

0.14999999999999990L  02 

-0.14999999999999990D  02 


0.53237900077244490D  01 
•0.16647580015448840D  02 
0.11999999999999990D  02 


N 


-0.65496972587504130D  00 

0.12395555984240360D  02 

•0.45626858219032430D  02 

0.42664577893729350D  02 


0.12290561919335860D  01 
-0.22746685660881600D  02 
0.76664299881176180D  02 
-0.57641919954206000D  02 

1  =  3 

-0.24952495419778680D  01 

0.42343845761147220D  02 

-0.96261459981441690D  02 

0.57641919954205890D  02 

1-4 

0.87783059330621940D  01 

-0.49135573227363170D  02 

0.82366875462154930D  02 

-0.42664577893729010D  02 


0.64457574561626680D  00 
■0.18660971287045070D  02 
0.11578146204084730D  03 
-0.239482314513344^0D  03 
0.15481333922979940D  03 

1  =  2 

■0.11196222405308090D  01 
0.32133165655799910L  02 

-0.19334776617762280D  03 
0.375637628382300801)  03 

-0.21703556145201720D  03 

1  -3 

0.187499999S9999980D  01 
-0.52499999999999990D  02 

0.28874999999999980D  03 
-0.47249999999999990D  03 

0.23625000000000010D  03 

1  =  4 

-0.37321558320700430D  01 
0.95791727360611900D  02 

-0.36864824974282230D  03 
0.49250461742576680D  03 

-0.21703556145201670D  03 

1  =  5 

0.13096091215873580D  02 

-0.11370836617381250D  03 

0.32621455387960370D  03 

-0.37977104240584490D  03 

0.15481333922979590D  03 
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ELEMENTS  OF  THE  EXPLICIT  INVERSE 
MATRIX  -WITH  DIVISION  BY  WEIGHTS 


N  =  6 

1  =  1 

-J.63670289144860900D  00 
0.26164498964931050D  02 

-0.24117639244238380D  03 
0.8234l6513ti9643150D  03 

-0.11597231598027260D  04 
0.57023452625729080D  03 

1  =  2 

0.10533062297523300D  01 
-0.43378702461460870D  02 

0.39436316819926010D  03 
-0.13115640254625750D  04 

0.17703674113121060D  04 
-0.81606492177331800D  03 

1  =  3 

-0.16209292953026820D  01 

0.65461713775432130D  02 

-0.57508919022082130D  03 

0.17945639576737070D  04 

-0.22082425392646350D  04 

0.92756392549416520D  03 

I  .  k 

0.26369381625442730D  01 

-0.10382467676719200D  03 

0.83478670215413420D  03 

-0.22372330555568110D  04 

0.24295770882061870D  04 

■0.92756392549416430D  03 

1  =  5 

-0. 521376395623408501)  01 
0.18819940606882500D  03 
■0.1w787736580489900D  04 
0.23907435979473020D  04 
■0.23099571975544600D  04 
0.81606492177330980D  03 


(n  =  6)  : 

1  =  6 

0.18277283982093760D  02 
.0.22634124784499930D  03 
0.97307045300343430D  03 
-0.188686913725838701)  04 
0.16914414714836860D  04 
■0.57023452625727700D  03 

N  =  7 

1  =  1 

C.63505671862000200D  00 

-0.34911537906429760D  02 

0.44427978535214690D  03 

■0.22113580870536630D  04 

0.50658074548755810D  04 

-0.53641853692382860D  04 

0.21240546327832390D  04 

1  =  2 

■0.10301560713220060D  01 

0.56505693626735990L  02 

■0.71390602118058990D  03 

0.35067953361694150D  04 

-0.78710473381123000D  04 

0.80943293742454660D  04 

-0.30785879539985880D  04 

I  =   3 

0.14678564240265540D  01 
-0.61203287902126410D  02 
0.10092970779827400D  04 
-0.48l31388752110o50D  04 
0.10337417338213540D  05 
-0.10039689543210950D  05 
0.35893498892627010D  04 


99 


ELEMENTS  OP  THE  EXPLICIT  INVERSE 
MATRIX  .WITH  DIVISION  BY  WEIGHTS 


.(».■■?) 


N  -.8 


•0.21874999999999980D  01 
0.11812500000000000D  03 
•0.1417500oOOOOOOOOOD  04 
O.63524999999999990D  04 
■0.12560624999999990D  05 
0.11261249999999980D  05 
-0.37537500000000020D  04 


0.35204555588710010L  01 
■0.18529521480577540D  03 
0.20377373884617530D  04 
-0.79266328307b75080D  04 
0.13979217961099200D  05 
■0.11496409792365190L  05 
0.35893498b92626960D  04 

I  «  6 

-0.69410653211762830D  01 
0.33499054543956470D  03 
■0.26553296088702320D  04 
0.86058593537967680D  04 
-0.13578219776863660D  05 
0.10377198349745970D  05 
■0.30785879539985680D  04 

1  =  7 

0.24321935531217310L  02 
•0.40620454164723460D  03 
0.24240160528103890D  04 
■0.68911106957305150D  04 
0.10105700100432590D  05 
■0.73801424274610320D  04 
0.21240546327832080D  04 


■0.632636290503882201)  00 
0.44904361134180810D  02 
■0.75130772461264850D  03 
0.50790321475309150D  04 
•0.16738927934245450D  05 
0.28517349619132120D  05 
■0.24099556251708940D  05 
0.79803684875954930D  04 

I  -  2 

0.10088738950655580D  01 
-0.71515d69608273520D  02 
0.11915715970486080D  04 
■0.79955698673191210L  04 
0.26057030694553270D  05 
■0.43705585221165210D  05 
0.36178290078771380D  05 
•0.11664144752831990D  05 


•0.14081395444254550D  01 
0.99539951323618810D  02 
-0.l6437571ti73o41720D  04 
0.10856214613981250D  05 
■0.34559345638578920D  05 
0.56179519216495980D  05 
-0.44749661063421930D  05 
0.13823425769235650D  05 

I  «  k 

0.19506271541102560D  01 
•0.13714860269817760D  03 
0.22260096024338430D  04 
-0.14261847226419680D  05 
0.43486763859943770D  05 
•O.D7093136369321480D  05 
0.50629390880600230D  05 
■0.14854809783306000D  05 
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ELEMENTS  OP  THE  EXPLICIT  IWV3RSE 

:  r.?AH     WITH     Dl/Iblo         ^    -       -vV:    ,:!io 


(N  ■  8) 

1  =  5 

-0.28270116134008710D  01 

0.1966206b3462510lOD  03 
-0.30804548507999980D  04 

0.18576680283564960D  05 
-0.524563971933704301)  05 

0.75267796535146170D  05 
-0.53354277602541830D  05 

0.14854809783305980D  05 

1  =  6 

0.45275221271127710D  01 
-0.30689695082148760D  03 

0.44110301906861850D  04 
-0.23240704879435370D  05 

0.50913236415C320400D  05 
-0.77973493S89913020D  05 

0.52014319321227460D  05 
-0.13823425769235560D  05 

1  =  7 

-0.89144666562789340IJ  01 

0.55415840227710460D  03 

-0.57814946771425290D  04 

0.25502564074450390L  05 

-0.58041610578821550D  05 

O.715828O455G008030D  05 

-0.45470723191052030D  05 

0.1l664144752o31820D  05 

1  =  8 

0.31230068535155660D  02 

-0.6756636l6u9512550D  03 

0.53201117676994600D  04 

■0.20618588633532940D  05 

0.43667373451622040D  05 

-0.5150775034&382530D  05 

0.31763023161456740D  05 

-0.79803684875953020D  04 


N  =  9 

1  =  1 

0.63094690159499350D  00 
-0.56144067147718930D  02 

0.1192^225317072080D  04 
-0.10446389121829060D  05 

0.4b224589234221030D  05 
-0.1121O523386107140D  06 

0.15126259454114080D  06 
-0.106139538456527&OI)  06 

0.30188269946895430D  05 

1  =  2 

-0.99434922630698750D  00 

0.8o40b27962o242270D  02 
-0.18725077842749480D  04 

0.1u335839433292440D  05 
-0.71786191232140630L  05 

0.17262944053493840D  06 
-0.22^9979883000793013  06 

0.158974559316330701}  06 
-0.443817C0059841140D  05 

1  =  3 

0.135612l70tt5072470D  01 
-0.12036985085193770D  03 

0.25359056o32479280D  04 
-0.21917518679328140D  05 

0.95003168642251510JD  05 
-0.22431467127417700D  06 

0.29207121661956960D  06 
-0.19644199b;2830700D  06 

0.53188530258027430D  05 

1  =  4 

-0.18076013195316310D  01 

0.159954125710033^013  03 

-0.33374743057706110D  04 

0.2b363935208562880I)  05 

-0.12000917873689330D  06 
0.27472299293684350U  06 

-0.34504174828970990D  06 
0.22333H7167123050D  06 

-0.58191387545852760D  05 
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(N  = 


ELEMENTS  OF  TEE  EXPLICIT  INVERSE 
MATRIX  -WITH  DIVISION.  BY  WEIGHTS 

9)  (N  «  9) 


I  =-5 

0.24609374999999980D  01 
-0.21656250000000000D  03 
0.44395312499999990D  04 
•0.36599062500000000D  05 
0.14850773437499990D  06 
•0.32376093749999990D  06 
0.38710546874999970D  06 
■0.23930156249999990D  06 
0.59825390625000010D  05 


■0.35423371992380150D  01 
0.30832612996597780D  03 
■0.61012780666898300D  04 
0.47404503239454910D  05 
-0.17882655111494290D  06 
0.36429059007332490D  06 
■0.41108239647497120D  06 
0.24219992849559070D  06 
■0.58191387545852670D  05 


0.56589921408547840D  01 

-0.48003575926795000D  03 

0.87216161808651660D  04 

-C.59461921498086270D  05 

0.20222o63H3608010D  06 

-0.38138919379832560D  06 

0.40625635414618650D  06 

-0.22906628353591180D  06 

0.53188530258027340D  05 

1  =  8 

-0.11134161372277090D  02 

0.8662o630780323l80D  03 

-0.11479412061649910D  05 

0.65739913426994570D  05 

■0.20121824117593400D  06 

0.35426794697368490D  06 

■0.35986367476130570D  06 

0.19607904116239290D  06 

■0.44381700059839920D  05 


0.39001694289452010D  02 

•0.10602788886162110D  04 

0.10622421003920810D  05 

■0.53508791074806880D  05 

0.15253238835015280D  06 

-0.257003142824798701)  06 

0.25355738385848900D  06 

■0.1353b66211l862080D  06 

0.30188269946892230D  05 


N 


10 


1  =  1 

•0.62972078889511680D  00 

0.6d631241594361980D  02 

-0.18007322514163060D  04 

0.19787480981637290D  05 

-0.11232455134290630D  06 

0.36273559299210630D  06 

-0.69145818794950720D  06 

0.76901847430852180L  06 

-0.46080521705457350D  06 

0.11482677561427350D  06 

I  =  2: 

0.98397308015031300D  00 

-0.107l8l875o6752850D  03 

0.28074636043399860D  04 

-0.30759373574475980D  05 

0.17385617851724700D  06 

-0.55819292467264940D  06 

0.1056149876743173013  07 

-0.11637950926080640D  07 

0.68955600212862640D  06 

-0.16952953248708180L  06 


1Q2 


iLEE.BNTS   0?    CE&   EXPLICIT   IN/^RSE 
MATRIX   WITH   DIVISION    BY    WS1  rHTS 


(N   =   10) 

1   =   3 

-0.132009-26200040000  01 

0.143638O6OO9058600O  03 
-0.37496162693665910D  04 

0.408401686761198600  05 
-0.22886025663o27430D  06 

0.726439281044946500  06 
-O.135466OO764180940O  07 

O.1467262o4050193300  07 
-0.05202751523799060D  06 

0.20479977286030410D  06 

1  =  4 

0.17134865755851150D  01 

-o.io6094715330446lOD  03 

0.48294021160310210D  04 
-0.5-06845G405179270O  05 

0.2o75415069b^2o020D  06 
-0.89548791099262U10D  06 

0.16320721886540260L  07 
-0.17216549991630^700  07 

0.97l838260913bb580D  06 
-0.22686995363410560D  06 

i  =  «> 

-0.-22391285740396210D  01 

0.24240619103871870D  03 
-0.62282226225405700O  04 

0.o60045o796l004blOD  05 
-0.35?7b41585769b580D  06 

0.107511267133415100  07 
-0.1o939045898428lOOD  07 

0.19286143921177980O  07 
-0.105l6'9213b889739OD  07 

0.23764031359797480E  06 


(N   =    10) 

1   =  6 

0. 30224413129666500D   01 

-0.32536632332934910D  03 

0.8212095  J859297140I)   04 

-0.84433163665506040D   05 

J. +36943bbbl372b230D  06 
-0.12605    7102770^0500   07 

0.212080260829882oOD  07 
-0.2C701235705269370D  0; 

0.108707066349202601;  07 
-0.237o4031359797450D  06 

1  =   7 

-0.43347755028589230D  01 

O.4ol52442191b95b900  03 
-J.11245195102679330D  05 

0.109076224196158500  06 
-0.526l96062d2311960D  06 

0.1423001308212d740O  07 
-0.22o6797605l63b040D  07 

0.21149872426799040O  07 
-0.10701713^179308300  07 

0.226089953634105600  06 

I  =  R 

0.D915269385903729OO  01 
-0.717^4Ob731312l890D  03 

0.1^0620363672135800  05 
-0.137129876124356000  06 

O.59ol79707423714900  06 
-0.150144898894764400  07 

0.22625302986972490O  07 
-0.202389454156694800  07 

0.99117-440504742600  06 
-0.20479977286030280O  06 
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ELEMENTS  OP  THE  EXPLICIT  INVERSE 
MATRIX  WITH  DIVISION  BY  WEIGHTS 


is  s  10) 

1  =  9 

■0.13600251472079610D  02 

0.12944478368158440D  04 

•0.21206736297286420D  05 

0.15243922236262340D  06 

■0.59948947934917960L  06 

0.14065755831522120D  07 

-0.20233^84408265790D  07 

0.17504102451139530D  07 

-0.83620979025508760D  0.6 

0.16952953248707710D  06 

I  m    10 

0.47636818940463100D  02 

-0.15888058270619700D  04 

0.19703318379258890D  05 

-0.12468504388204750D  06 

0.45693572875155430D  06 

-0.10264559981162400D  07 

0.14345742062O0863OD  07 

-0.12163406599856000D  07 

0.57263576347380240D  06 

-0.11482677561425580D  06 

N  -  .11 


(N-11) 


■0.97629228872923490D  00 

0.12783586044824230D  03 

-0.40528073056744170D  04 

0.543407541235H910D  05 

-0.38217790816757760D  06 

0.15651228711533680D  07 

•0.39233b53943289440D  07 

0.60947715541134540D  07 

-0.57158374375782260D  07 

0.29608725820561420D  07 

-0.649ol738642927510D  06 

1  =  3 

0.12940094492703080D  01 
-0.16931341404173450D  03 

0.53555796829995070D  04 
-0.71527329877500880D  05 

0.50020326653883290D  06 
-0.2033H29052140&90D  07 

0.50485603123939950D  07 
-0.77536105791211750D  07 

0.71747065316893470D  07 
-0.36600707958978520D  07 

0.78967223585514890D  06 


0.62880260230913370D  00 
■0. 8236622. 639663810D  02 
0.26142904642913460D  04 
■0.35122822330837650D  05 
0.24773476663498850D  06 
•0.10184725075194570D  07 
0.25656309060809040D  07 
■0.40097677383656480D  07 
0.37879189896732840D  07 
•0.19791458507143100D  07 
0.43874883894089130D  06 


■0.16477565202582030D  01 
0.21533446996204010D  03 
■0.67853720573359590D  04 
0.90030823198171330D  05 
■0.62369982976702370D  06 
0.25040530553528440D  07 
-0.61246301487788360D  07 
0.92410882893189820D  07 
■0.83826594326138520D  07 
0.41852673378896550D  07 
-0.88288361424753310D  06 
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ELEMENT.     OF    THE   EXPLICIT    INVERSE 
MATRIX    »VITH   DIVISION    BY    HEIGHTS 


(N  =   11) 

1  =  5 

0. 209253200, 9^u4940D  01 

-0.27291771296259 iOOD  03 

0.854701597110103400  04 

-0.112<:12761080o8450C  06 

0.76580444244534^30D  06 

-0. 3-1633764549051100  07 

o.7^126476877351980D  07 

-0.10613016080694200D  08 

0.937660396353-'6540D  j7 

-0.4558773066o363300D  07 

0.9370111O644ob6970D  06 

1  =  6 

-0.2707031250D000030D  01 

0.35l9140o249999990D  03 
-0.10909335937500C10D  05 

0.14076562500QOOOOOL  06 
-0.93785097656^50140D  06 

O.358713O421&75002OD  07 
-0.83037642187500140L  07 

0.118214971875o0010O  08 
-0.101164735546O7500I)  08 

0.47740661716750060D  07 
-0.95481323437500150O  06 

1=7 

0.3636o46o696400980D  01 

-0.47010605445151210D  03 

0.14314916797832990L  05 

-0.179234601267.78200  u6 

O.1147471006292O09OO  07 

-O.41974249472410480O  07 

0.9301d40o412828020D  07 

•0.12725317995043420D  08 

0.10513146152103770D  00 

-0.46113379976303560D  07 

0.93701110644666O20T;  06 


(N  -  11) 

I  =  ft 

-0. 5204^9140348259700   01 

0.66^41^173434762700  03 
-0.195620738195873COJD  05 

0.231256345102216901)  06 
-0.13O20313593756O70O  07 

0.47527C82042786310D  07 
-0.999457884598442300  07 

0.13096596717268040D  08 
-0.104450160327458600  08 

u.464356o8045856490L  07 
-0.88288361424752910D  u6 

1-9 

0.0^  :oo41;111262ol90D  01 
-0.10  $36c  ';9056o76870D  04 

0.27931^099685404300  05 
-0.29127968457733320D  06 

O.1578379621O432000D  u7 
-0.50444757345765650D  07 

0.10O50291020008940D  00 
-0.1264^1613246886000  00 

0.97693199820902740O  07 
-0.423665156265359700  07 

0.789672^3585514^00D  06 

I   =   10 

-0.l6312795060o2b4oOO  02 

0.16644O68765925460O  04 
-O.3b96524o4451012300   05 

0.3^5U164^79o71olOO  06 
-O.159036o5l491343500   07 

0.4755/994O39214340O  07 
-0.90517670251564250O  07 

0.1101O601364003760O  08 
-O.03097o658039000200  07 

0.35373012622364760D  07 
-0.64-J01738642925160O  06 
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ELEMENTS  OF  THE  EXPLICIT  INVERSE 
MATRIX  WITH  DIVISION  BY  WEIGHTS 


(N  «  11) 
I  *  11 

0.57135446067443050D  02 

-0.2<±931949503715550D  04 

0.34449271098960500D  05 

-0.2669H0824436b290D  06 

0.12171732421673420D  07 

-0.34859840680^22630D  07 

0.64479931659575390D  07 

-0.76941942262113810D  07 

0.57193040855838490D  07 

-0.24083425386942900D  07 

0.43o748d3894083490D  06 

N  «  12 

I  -  1 

-0.62809712046008180D  00 

0.97349208725481830D  02 
-0.36740962081804510D  04 

0. 5916423620406 3700D  05 
-0.50613790704194580D  06 

0.25673908975000990D  07 
-0.8l82553'ol31982310D  07 

0.16806961478159360L  Oo 
-0.22201731817903030D  08 

0.18207319461823990D  08 
-0.84295810775423920D  07 

Q,1682bl321467249$OD  07 

1  =  2 

0.97044216519729790D  00 

-0.15036966861871660D  03 

0.56705340420321710D  04 

-0.91l85184277295b30D  05 

0.77849448d935539dOD  06 

-0.39383401739560010D  07 

0.12509469255818680D  08 

-0.25538333143220000D  08 

0.3363522C804863140D  08 

-0.2742431589034b060D  08 

0.12611867398301090D  08 

-0.24984179627599490D  07 


(N  «  12) 

'  I  ■-  3 

■0.12744705976447860D  01 

0.19737725411053610D  03 
-0.74314329933348770D  04 

0.11917694554194690D  06 
-0.10134999401397680D  07 

0.51007872028815800D  07 
-0.16097249175721700D  08 

0.32670693022309940D  08 
-0.42551737919709900D  08 

0.34329427598729830D  08 
-0.15600219960168470D  08 

0.30498673596812220D  07 


0.15997929592334350D  01 

■0.24755198318340960D  03 

0.92964442261903590D  04 

■0.14842587620106170D  06 

0.12542388006541950D  07 

-0.62600943174096300D  07 

0.19553828680241450D  08 

-0.39205841039798720D  08 

0.50356479548145220D  08 

-0.40001461744518260D  08 

0.17875659251661460D  08 

-0.34334399481679000D  07 

1  =  5 

-0.19896588o93697860D  01 

0.30747756727537110D  03 
-0.11500307185105370D  05 

0.18234946701209440D  06 
-0.15258580461689860D  07 

0.75200497695791070D  07 
-0.23133310008012750D  08 

0.45577110459247d80D  08 
-0.57424066378613350D  08 

0.44696839423953880D  08 
-0.19561620276687970D  08 

0.36797047140185570D  07 
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li  3    JKTb   OF    THE    EXPLIC  I  C    [tfVBSSE 
LA.1; IX    tflTH   DIVISION  ,      I  I     PS 


(N  =   12) 


0.24981120367529350D  01 

-J.3b526531093196790D  03 

,;.1431b8o68590b215oO  05 

-0.2^46l049t)37321410D  06 

0.1b51324413«972jdOD  07 

-O.b95242o5939ib02300  07 

0.2693o33iSl425B800D  08 

-0.510124627094654100  08 

o.b36782l2531130^o  Q'd 

.0.48357186859487330D  08 

0.20665514635^393000  Ob 

-0.380063598dl2573oOD  07 

1  =   7 

-0.32133818890451120D  01 
0.493940742022827900  03 

-o. 181702629695564100  05 

0.2b01164157212b410D  06 
-O.22544988l^6o099100  07 

0.10593935323414500L  08 
-0.30O98406193109330D  08 

0.575936093587351200  08 
-0.686932436652302100  06 

0.507370196560090500  08 
-0.211414dl234343w90D  08 

0.3o00635968l2573bOD  07 

1  =  8 

0.430505169011563800  01 
-O.657968ll3o49961100  03 

0.2377451766oo51740O  05 
-O.355719911o2c5bb500  06 

0.275302519211o250D  07 
-0.1236682773^7401200  08 

0.34650435989768490D  08 
-0.6217892111d244410D  08 

0.717256^37990730800  06 
-0.514643959260942000  08 

0.209151315775160100  08 
-0.36797047140165370D  07 


(N  =   12) 

1   =   9 
■0.6153  151  3164492020D  01 

0.9301 jc 44277o089  30D     03 

•0. 32449465759824320D  05 

0.45  0097275 3631060D  06 

•0.33'--05904613610080D   07 

0.14064233ol9780170D  08 

■J. 373  .92309944750200  08 

0.b436270015332.790O  06 

-0.71769601275455330D  Oo 

0.50O64066377137580O   08 

■0. 193921 O0178185280D  06 

0.34 334 3994blb7b820D   07 

I   =   10 

0.98C3194840362150oL  0. 
-0.14440o92792927150D  04 

0.46  ^26001574812040D  05 
-0.57o63539694709920D  06 

0.3^0221433097761400  07 
-0.149956097635312300  08 

0.37613307059219770D  Ob 
-O.6254601cu,934b98200  06 

O.D7631326608678030D  06 
-0.460699327795117900  08 

0.1794o3209963247700  08 
-0.30490673596M1O90O   07 

I  =  11 

-0.1^271027^76777.:-50O  02 

0.26O44176581899O1OO  04 

-0.614215396866676100  05 

0.6479296405&3506600  06 

-0. 38469474672530560O  07 

0.142064721bl404030D  08 

-0.34^4216O151b55240O  06 

O.54b3577oo714935000  Ob 

-0.57888553140108210U  08 

0.38718629859131920D  Oc 

-0.14870730192057450D  08 

0.24964179627597960D  07 
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ELEMENTS  OP  THE  EXPLICIT  INVERSE 
MATRIX  WITH  DIVISION  BY  WEIGHTS 


I  «  12 

0.67497577870248550D  02 
-0.320859l82bb8901bOD  04 

0.573750304b6011780D  05 
-C.5334051611/338810D  06 

0.29535-3612699131^  07 
-0.10448804284924580D  08 

0.24480199520638750D  03 
-0.38435^39097431290D  08 

0.39997175270061230D  08 
-0.26466235493383090D  08 

0.10081364283853400D  08 
-0.16828132146722190D  07 

N  =.13 

1  =  1 

0.62754331553197030D  00 
-0.11358033610843180D  08 

0.50.-50939417970050IJ  04 
-0.95442184513324400D  05 

0.97168744238521400D  06 
-0.59404968247014980D  07 

0.23229631009714890D  08 
-0.60063715334341270D  08 

0.10380582544494940D  09 
-0.11855578799997500D  09 

0.85833H7642566710D  08 
-0.35664862301875860D  08 

0.64752096878575660D  07 

1  =  2 

-0.96588072089159760D  00 

0.17478290345705200D  03 

-0.7728269474^5b3250D  04 

0.14663509169035800D  06 

-0.14907030198521330D  07 

0.90960216544733470D  07 

-0.354&3079564821880D  08 

0.91477792834819540D  08 

-0.15754736564692130D  09 

0.17920322086746790D  09 

-0.12913513896815610D  09 

0.5337199515cH17440D  08 

■0.96318464320360660D  07 


(N  =  13) 

I  -3. 

0.125942b6332923030D  01 

-0.2278178691145bl30D  03 

0.10061811146213320D  05 

•0.19054009023605280L  06 

0.19316439162564420D  07 

-0.1174334852«528470D  08 

0.45600386924790600D  08 

-0.11691192072352450D  09 

0.20004546717896550D  09 

-0.225b4442129174520D  09 

0.16137022920607280D  09 

-0.66066587924320450D  08 

0.117992b7514586200D  08 

1  =  4 

-0.15635930378826840D  01 

0.28266984060442140D  03 

-0.12461600812048860D  05 

0.23524972654785430D  06 

-0.23742867586914440D  07 

0.14350435738888250D  Ob 

-0.55323068110608670D  08 

0.14062525506500070D  09 

-0.23o24303520694750D  09 

0.265975476323290401)  09 

-0.18771579270607220D  09 

0.75836132834923010D  08 

-0.13354193591b6l380D  08 


0.1914224o647590760D  01 
-0.34574583103880170D  03 

0.15200112938246720D  05 
-0.2b559473259563270D  06 

0.28631024340339590D  07 
-0.17154969863530380D  08 

0.65436802612953100D  08 
-0.1642b685625950120C  09 

0.2744bl06194385310D  09 
-0.30181562617390290D  09 

0.209o0921505546b20D  09 
-0.832bl068929590510D  08 

0.14419082659029250D  08 
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3Li&   i  .        Oj*    CHE   EXPLICIT    LtSTVii  i:  J 

...;.MI)t   WITH  DIVISION    BY    .'/EIGHTS 


(N_=    13) 

1  =  6 

-0.23511020637611940D  01 

0.4240790o93(J1'j4dOOj  03 
-0.l85662232b504oOOOO  05 

0.34637477026633250D  06 
-0.34376691560332u90D  07 

O.203376910794b9o200  08 
-0.76406714822997720D  08 

u.ld85bl468469960oO!;  09 
-0.309244140bl3ol7  30D  09 

^.333533203110b4'j30D  09 
-0.22715895046998290D  09 

0.d853H9829b31113OI>  Ob 
-0.15044119928021750D  Ob 

1  =  7 

0.2^326171875000020D  01 

-0.527b7109375000080D  03 

0.22962392578125020D  05 

-0.42376318359375070D  06 

0.41429083007olL5oOO  07 
-O.^405220004oo749b0.i.;  08 

0.88433892421875160D  08 
-0.21324162234375030D  09 

0.34154465071289120L  09 
-0.3599H5441796bbl^O  09 

0.23973768959765660D  09 
-0.91502934960937660L  08 

0.152504891601562oOD  08 

1  =  8 

•0.37592935460216400JD  01 

0.67442121115^569700  03 

-O.2903582744^o99610j  05 

0.5266^500414.  8^910D  06 

•0.50^851593966171001:  07 

0.28384930929296120D  08 

■0.101256110701618801)  09 

0.^3o91794u05c^3040L  09 

■0.36b789729244o4/20O  09 

0.37854677945665700D  09 

-0.24622768445999690O  09 

O.9199824uo399498900  08 

•0.150441199200217200   Ob 


(N  =  13) 

E  =  9 

0.5-275686326371340D  01 

-O.89b7o5387598239000  03 

0.3792 3661623981460B  05 

-0.667ob80552b9478lOO  06 

0.61340924932227670'D  07 

-0.33194249772054250D  08 

0.1136996 J774122740L  09 

-0 . : ■  564656 2 1 67570360D  09 

0.3b662464671274010D  09 

-0.386015918239934301)   09 

0.24517691232592170D  09 

-0.89747922978760040D  08 

0.14419082659029160D  08 

I  =  10 

-0.718009578502850001)  01 

0.12666259857963150D  04 

-0.51721949968460710D  05 

0.b6073246'5039bl2O0I)  Ob 

-0.74051109220683070D  07 

0.37765288364222390E  Ob 

-0.1-  )iO:j23471440b30D  09 

0.2666944929104176OD  09 

-0.38903833027967030D  09 

O.37cll77  3502ol67300  09 

-0.234b951085b4767900  09 

0.84414190267412770D  08 

-0.1335  193591861230D  Ob 

I  =  11 

0.11435022541389190D  02 

-0.19659144^59bl4bOO!.    04 

0.738373370890617 300  05 

-0.108718740227405901)  07 

0.85006375595485470D  07 

-0.40411727534986120D  Ob 

0.12^06633008148290D  09 

-0.260526011460648701)   09 

0.36875640203382960D  09 

-0.35003438bl4031490O  09 

0.21338941800123120D  09 

-0.75524o22250711790D  Ob 

0.117992675145b5840D  08 
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ELEMENTS  OF  THE  EXPLICIT  INVERSE 
MATRIX  WITH  DIVISION  BY  WEIGHTS 


(N  «  13) 

I  M   12 

•0.22477370421262190D  02 

0.35453247719126610D  04 

■0.98039763299629560D  05 

0.122041596b8l98080D  07 

-0.d6285698775474b90D  07 

0.38431084088930520D  08 

•0.11373773829282280D  09 

O.22946582918523830D  09 

-0.31718441412569740L  09 

0.29569465014905760D  09 

■0.17774505673993810D  09 

0.6221Ol62O2bO1317OD  08 

•0.96318464320357350D  07 

1-13 

0.7872321576173544OD  02 

•0.43733370845259480D  04 

0.9174924<^665078710D  05 

■0.1U069243767794210D  07 

0.66410657217231400D  07 

■0.28341539376072510D  08 

0.81542617182726970D  08 

•C.16131015049826300D  09 

0.21982054304058640D  09 

-0.20^75409315112^90D  09 

0.12088347172049520D  09 

■0.42037653952402330D  08 

0.64752096o7b55o300D  07 


N  »  14 

I  «  1 

■0.62710058676849130D  00 

0.13U5968903823980D  C3 

■0.67159735263579190D  04 

0.14d45758311417680D  06 

■0.17713639267792060D  07 

0.12bl4838745636220D  08 

■0.60088842697869890D  08 

0.18975509790027530D  09 

■0.41111999022164400D  09 

0.61198659250115030D  09 

-0.61492202573029410D  09 

0.39829832256181480D  09 

■0.15007902306178370D  C9 

0.249O4612699678530D  08 


0.962253563152O779OD  00 

■0.20107526446341820D  03 

0.10299-71598740270D  05 

-0.22749551893141230D  06 

0.2711502407462771CD  07 

■0.195b8419447538950D  08 

0.91687380045722230D  08 

■0.28891943525264980D  09 

0.62438119257867880D  09 

-0.92670839285596560D  09 

0.92801382186042630D  09 

-0.59879739664643370L  09 

0.22465816840735830D  09 

-0.37221050666135830D  08 
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ELEMENTS   OF   THE   EXPLICIT   INVERSE 
MATRIX   WITH  DIVISION    BY    WEIGHTS 


(N  «   14) 
1  =  3 

•0.124758652l27b4l20D  01 

0.260627b07177b7190D  03 

•O.1333837OO6625103OD  05 

0.294202«443612.,140D  06 

-0.3499406b475626250D  07 

0.25212540726blO«10D  08 

-0.1176l7540ol302880D  09 

0.36913741442143300D  09 

■0.79397229603140680D  09 

0.117201118S043305OD  10 

•0.11664350161564310D   10 

0.74745627373425O50D  09 

-0.27839900529629120D  09 

0.45724736356779090D  08 

1  =  4 

0.1535528b284604530I)  01 

-0.32064094O43427940D  03 

0.163879094o2115830D  05 

■0.36064900873599210D  06 

0.42759286645134930D  07 

•0.30677579805035340D  08 

0.14236676990651250D  09 

-0.44403634576476030D  09 

0.94819133467043650D  09 

•0.13882270207062360D  10 

0.13690686886524820D   10 

■0.86858664432220080D  09 

0.31993822789747540D  09 

■0.51968787273840710D  08 


( N  =  14 ) 

1  =  5 

■0.18570160912215590D  01 

0.3875222741^995950D  03 

•0.19767119057bb677013  05 

0.43355769719^399801)  06 

-0.5115b9221778lb850D  07 

0.36477408002075480D  08 

■0.16800267212292770D  09 

0.51933559363838030D  09 

-0.10977565123736930L   10 

0.15b91593037308090"D  10 

-0.15481701995747850D   10 

0.96954303861109500D  09 

■0.35232510694594320D  09 

0.5b44086b8l4U93620D  08 

1  =  6 

0.^2429767477604370D  01 
•0.46762438068047850D  03 

•J.23784307536616040D  05 
■0.5191205259139764CD  06 

0.6u83296839682C860D  07 
-0.42991511298815780D  08 

0.1958916643129211013  09 
■0.59bl06097577o9330T;  09 

0.12470305777461690D  10 
■0.17788305398299360D  10 

0.170b439269b70oObOD  10 
-0.105196l58049340oOD   10 

0.3762b797956669640D  09 
■0.593472&1881679930D  Ob 
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ELEMENTS  OP  THE  EXPLICIT  INVERSE 
MATRIX  WITH  DIVISION  BY  WEIGHTS 


(N  «  1*0 

1  =  7 

•0.273475134358l8^D  01 

0.56936165211156270D  03 

■0.2883580699364rf890L  05 

0.62485803994337690D  06 

■0.724886077593ol210D  07 

0.50579039189955630D  08 

-0.22702150561159Q70D  09 

0.681599522641658001;  09 

•0.13958347470868530D  10 

0.19546601596146270D  10 

•0.18408185547487770D  10 

0.11145041508177720D  10 

•0.39179754849029220L  09 

0.607S1758986943110D  08 


8 


0.33973558713789260D  01 

■0.70582687361741650D  03 

0.35516457151659940D  05 

-0.76127530737318370D  06 

0.86993044124879000D  07 

■0.59581548577651650L  08 

0.26184043920873030D  09 

•0.76867050936797620D  09 

0.15387939241834260D  10 

-0.21076999149043930D  10 

0.19436719741531690D  10 

■0.11539107699158290D  10 

0.39836531833996820D  09 

-0.60781758986943040D  08 


(N  »  14) 

1-9 

•0.43454135994380020D  01 

0.8997728889565H30D  03 

-0.44808975591674030D  05 

0.94395747735251700D  06 

•O.1O5375364O8041470D  08 

0.70205644801013560D  08 

■0.29957197201463970D  09 

0.85419165997620520D  09 

-0.166371407416518801)  10 

0.22222615089891290D  10 

-0.20034483623153590D  10 

0.11655922524640850D  10 

-0.39522642489514270D  09 

0.59347261881679870D  08 

I  -  10 

0.58047100121171420D  01 
-0.11950398764140700D  04 

0.5b456906l95874110D  05 
-0.119559877^1625510D  07 

0.128473420064b0540D  08 
-0.82123426856577030D  08 

0.33681964482398110D  09 
-0.92684340544247790D  09 

0.17500898366654430D  10 
-0.22760221055625940D  10 

0.20054346475457740D  10 
-0.11440295227590730L  10 

0.3&1406lb763727l60D  09 
-0.5b4408688l4093220D  08 
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3LEMi*NTS  OP    THE   EXPLICIT    INVERSE 
MATRIX   WITH  DIVISION    BY    WEIGHTS 


(N  =   Ik ) 
I   =   11 

-0.b2b534743267l86bOD  01 

0.16d6931b54741400uD  04 
-0.796b67b0156054520I)  05 

0.15411315559^972201)  07 
-0.15521oOd435712500D  08 

0.9360295193111o650D  08 
-0.3b572146347742750L  09 

0.9673908689^435690D  09 
-0.17689635b76o97900D  10 

0.224107b3352539520D  10 
-0.19325345179767670D  10 

0.1C82b933169120550D  10 
-0.35565600666._44880D  09 

0.51966787273839930D  08 

I  =  12 

0.13192162617824440D  02 
-0.26176698bb9117130D  04 

0.11376129b41733790L  06 
-0.1^4b6373264239120L  07 

0.178546024869962^0D  08 
-0.10045213191534110L  09 

0.37287171793200110D  09 
-0.94b96563bb329b660D  09 

0.16845659139o204bOD  10 
-0.20b5l6l4127761090D  10 

0.176512424340394oOD  10 
-0.97439764600749770D  09 

0.31612^56734182530D  09 
-0.45724736356777360D  08 


(N_=J,4) 

I  =  13 

-0.25929439201462620D  0^ 

0.4/205431754915600D  04 

-0.151^26455000o0860D  06 

0.219179309431b6830D  07 

-0.181703146097125201)  08 

0.958199654365879901)  Ob 

-0.340250^6907390060D   09 

0.b3b907b71971312bOD  09 

-0.14P46416517418060D   10 

0.17686bll664737440D   10 

■0.147653b9l6b7943501j  10 

0.806141327716665O0D  09 

-0.2592154902523bb40L  09 

0.372^1050666132720D  08 

I  ■  14 

0.908l2360ob3i00790D  02 

-0.582b966390b486240D  04 

0.14173147839979980D  Oo 

-0.1bll5731489670740D  07 

0.14012562197212460L  08 

■0.70bl4860274995020L  Ob 

0.24448o590251506bOD  09 

-0.59112499164761420D  09 

0.101U5715360278980D  10 

-0.12157870827753120D  10 

0.10067432324796280D  10 

-0.54ol49b3639513160D  09 

0.17472094203397150D  09 

■0.249b461269966912013  08 
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ELEMENTS  OF  THE  EXPLICIT  INVERSE 
MATRIX  WIT*  DIVISION  BY  WEIGHTS 


N  «  15 
I  *  1; 

0.62674114724037670E  00 

■0.14978732641334850D  03 

0.879917021581786601)  04 

■0.22384410147813O90D  06 

0.30907278460131640D  07 

-0.26070735053359J80D  08 

0.14398922790202560D  09 

-0.54295844801453790D  j9 

0.14312370236968300D  10 

■0.26615478102941000D  10 

0.34745206974963210D  10 

-0.3H423495852369bOD  10 

0.1824261381267898GD  10 

■0.62b70665172053600D  09 

0.9663484325400bl50D  08 


0.l2380d714595348lOD  01 

-0.295801295523878201)  03 

0.17359797991491830D  05 

•0.44089167033424270D  06 

0.60732983ol4250210D  07 

■0.51071425638271780D  08 

0.2809879664i)939410D  09 

•0.1054674685H27840D  10 

0.27650b29106754510D  10 

-0.5109938155C103200D  10 

0.66235585027415350D  10 

-0.58895796177453840D  10 

0.34195573387889630D  10 

-0.1167O455028865680D  10 

0.17747321082681620D  09 


1  =  2 

-0.95932021867915070D  00 

0.229246611C7392350D  03 

■0.13462431736329^601   05 

0.34227829881890100D  06 

-0.472^1560871994830D  07 

0.39789495277815740D  08 

-0.21946537164326570D  09 

0.82623583095249640D  09 

-0.2173832422737238OD  10 

0.40336319906660950D  10 

-0.52525395858817300D  10 

0.469458655202855701)  10 

-0.27413061611607320D  10 

0.94143312509001210D  09 

-0.14414037028721530D  09 


■0.15132903440408250D  01 

0.361434582777l8l^OD  03 

-0.21190591940268440D  05 

0.53727739053594980D  06 

-0.73832858170019010D  07 

0.61893331012305830D  08 

-0.339211944026647001)  09 

0.12673395948070550D  10 

-0.33048142182154370D  10 

0.60700974692866130D  10 

-0.78144193438138940D  10 

0.68962211361328180D  10 

-0.39712890466219770D  10 

0.13434429202464070D  10 

-0.20239306917933810D  09 


114 


ELEMENTS   OP    THE   EXPLICIT    INVERSE 
MATHIX   WITH  DIVISION   BY    fffilG  ITS 


(N  *  15) 

1  =   5 

0.18124564883o14930D 

■0.432u621065176l890D 

0.2433103129b391310D 

-0.64067995696336^903 

0.877358790299329oOIj 

-u. 732156450625498703 

0.399036610113495103 

-0.14bl06bb7607554903 

0.3b330647941031090D 

-0.69d^95564132o2740D 

0.o903794ct)16235820D 

-0.777914853673266103 

0.44322964323493450L 

-0.14o27963179240750D 

0.22uc 344992o -79640D 

1  =  6 

-0. 216071642  3482^8701) 

0.51547223690100050D 
-0.30115551313185460L 

0.759028066159013603 
-u.lu34291371825l820D 

0.85762161473356650D 
-0.46 3800455 12426430L 

^.170592788286904703 
-0.43701651364934^903 

0.787056840743666903 
-0.991694025 3391^7003 

0.655600556664314703 
-0.481351289631718703 

0.158940354717114903 
-0.233635341347084203 


(N  -  15) 


1  =  7 

01 

O.2^bo95o50657809103 

01 

03 

-0.61703749015-204403 

03 

05 

0.359428395360953103 

05 

06 

-0.901380357017199903 

06 

07 

0.121966567554737003 

Ob 

Ob 

-0.100230128929750703 

09 

09 

v.. 536252329650570103 

09 

10 

-^.194b37996342654303 

10 

10 

0.49^450c62010870303 

10 

10 

-O.b742ob55747 3741303 

10 

10 

0.10856370799125670D 

11 

10 

-0.922764710195435503 

10 

10 

0.5113565^4467012403 

10 

10 

-0.166409954235023503 

10 

09 

0.241194716946680403 
1  =  8 

09 

01 

-0.31420698437500O603 

01 

03 

0.747817382812501603 

03 

05 

-0.433734082031250703 

05 

06 

0.107964830078125003 

07 

08 

-0.1446428  38183594003 

08 

08 

0.117362460056593803 

09 

09 

-O.6I06 39408105470003 

09 

10 

0.221100969375000303 

10 

10 

-0.549161007863770403 

10 

10 

0.957839651323243903 

10 

10 

-0.11686231979296680D 

11 

10 

0.976467034511720803 

10 

10 

-0.532351003754883703 

10 

10 

0.1705876l4462b90803 

10 

09 

-0.243696592089844203 

09 
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ELEMENTS  OF  THE  EXPLICIT  INVERSE 
MATRIX  WITH  DIVISION  BY   WEIGHTS 


(N  «  15) 


0.3^931009622298320D  01 

■0.92459692959775470D  03 

0.53278534088889670D  05 

-0.13120350817553940D  07 

0.17312394046906410D  08 

■0.13791130423589840D  09 

0.712O3877735789710D  09 

-0.24897282238044170D  10 

0.60497726835068710D  10 

-0.10329636183294670D  li 

0.12352321377315810D  11 

-0.1013048»499360400D  11 

0.54290104362649800D  10 

-0.17126264949032910D  10 

0.24119471694668040D  09 

I  -  10 

-0.49721421265370650D  01 

0.11769003011083430D  04 

■0.67H7059857607340D  05 

0.16244931453902490D  07 

-0.20943608396814680D  08 

0.16235745200514600D  09 

■0.81442361197469660D  09 

U.27681423968320320D  10 

■0.65498909091008840D  10 

0.10915782300378510D  11 

■0.12772292374733760D  11 

0.102739367601519C30D  11 

-0.54120828456768640D  10 

0.16814912316880170D  10 

■0.23363534134708260D  09 


(N  ■  15) 
I  4j.11 

^.66366757805783670D  01 
-0.156l86398w7458830D  04 

0.874907b5271139800D  05 
-0.20562089908785320D  07 

0.25527491884J02960D  08 
-0.1899921bl92076590D  09 

0.91677654852783980D  09 
-0.30096847588033260D  10 

0.6909589944106407010  10 
-U.11220336871860480D  11 

0.1284031212659236CD  11 
-0.10134u53583039l80D  11 

0.5251883732250801010  10 
-0.160888667175506501)  10 

0.22083449926279530D  09 

,  I  »  12 

-0.94692047715028570D  01 

0.22038877080534270L  04 
-0.1192^5b35506l4790D  06 

0.26504614225451120D  07 
-0.308fo3745176710530D  08 

0.216882744402032^UD  09 
-0.99775745931098660D  09 

0.315O8987956921740B  10 
-0.70097464706489320D  10 

0.11094760259234770D  11 
-C.12431850981447720D  11 

0.96417768253899200D  10 
-0.49243003787383290D  10 

0.14900600482642920D  10 
-0.20239306917933360D  09 


116 


ELEMENTS  OF  THE  EX  FLIC  IT  INVERSE 
Iv:/'  THIX  WITH  DIVISION  BY  HEIGHTS 


(N  =  15) 


(N  =  15) 


I  =  13 

0.13-74711637672710D  02 

-0.34193093774842380D  04 

0.170214738094u3o90L  06 

-0.3354^232418030070^  07 

0.35561798593i52b060D  u8 

-0.23330092608529340D  09 

0.10202340^66^0,:  310D  10 

-0.31^1350410^6745900  10 

0.67J05-7e094155060D  10 

-0.103652159413701^00  11 

0.11404637279697670D  11 

-0.8715807963530524^D  10 

O.4398027986503o2301)  10 

-0.131757944860C78SOD  10 

0.177473210b2680b70D  09 


I  -  ik 

0.10376501328106770D 

-0.76202104731309220D 

0.21252068208342310D 

■0.31278310698445420D 

0.2c0133584l6827970D 

-0. 16516,;77951213570L 

O.672032557768883OOD 

-0.19413456172002270D 

0.40403272170379310D 

■0.607594053914125800 

0.654056302260^58600 

-0. 491286  i;72o9458650D 

0.24448456450144930D 

-0.72418115  3835-2920-0 

0.96634843253961540D 


03 
04 
06 
07 
08 
09 
09 
10 
10 
10 
10 
10 
10 
09 
08 


1  =  1^ 

-0.2962d043653750330D  02 

0.6l66058l9bb535o80D  04 

-0.2_64&7421u5995170D  06 

0.37788653936583870L  07 

-j. 36266797135704260O  08 

0.223C9270240057440O  u9 

■0.9335-f3404^0152000D   09 

0.2749847&339129390L  10 

-0.580429454293298800  10 

0.882119583758752300  10 

•0.95729310319341260O   10 

0.72363984094249240O  10 

-0.36194492311267030D  10 

0.10765320589308O60D  10 

-0.14414u37028719950O  09 


117 


INITIAL  DISTRIBUTION  LIST 

No.  Copies 

1.  Defense  Documentation  Center  20 
Cameron  Station 

Alexandria,  Virginia  22314 

2.  Library  2 
Naval  Postgraduate  School 

Monterey,  California  93940 

3.  Commander,  Naval  Air  Systems  Command  1 
Navy  Department 

Washington,  D.  C.   20360 

4.  Professor  D.  J.  Collins  (Thesis  Advisor)  5 
Department  of  Aeronautics 

Naval  Postgraduate  School 
Monterey,  California  93940 

5.  Chairman,  Department  of  Aeronautics  1 
Naval  Postgraduate  School 

Monterey,  California  93940 

6.  Professor  A.  E,  Fuhs  1 
Department  of  Aeronautics 

Naval  Postgraduate  School 
Monterey,  California  93940 

7.  Ensign  Terrill  Jo  Wendt,  USN  3 
1932  Hempstead  Road 

Rock  Hill,  South  Carolina  29730 

8.  Dr.  E.  S.  Lamar  (Code  03C)  1 
Chief  Scientist,  Naval  Air  Systems  Command 

Navy  Department 
Washington,  D.  C.   20360 

9.  Dr.  R.  S.  Burington  1 
Chief  Mathematician,  Naval  Air  Systems  Command 

Navy  Department 
Washington,  D.  C.   20360 

10.  Commander,  Naval  Ship  Systems  Command  1 
Navy  Department 

Washington,  D.  C.   20360 

11.  Office  of  Naval  Research  1 
Navy  Department 

Washington,  D.  C.   20360 


118 


12.  Commander,  Naval  Ordnance  Systems  Command 
Navy  Department 

Washington,  D.  C.   20360 

13.  Commander,  Naval  Ordnance  Systems  Command 
(Attention:   P.  D.  Fisher,  Code  03211) 
Navy  Department 

Washington,  D.  C.   20360 

14.  Commander,  Naval  Ordnance  Systems  Command 
(Attention:   B.  E.  Drimmer,  Code  033) 
Navy  Department 

Washington,  D.  C.   20360 

15.  Superintendent 

U.  S.  Naval  Academy 
Annapolis,  Maryland   21402 

16.  Head,  Engineering  Department 
U.  S.  Naval  Academy 
Annapolis,  Maryland  21402 


119 


UNCLASSIFIED 


Security  Classification 


DOCUMENT  CONTROL  DATA  -  R&D 

(Soourity  cloooitlcotlon  o/  tttto,  body  ot  mbttrmct  ond  Indoxtng  annotation  mutt  bo  ontorod  wftan  tho  ormrall  roport  to  elmootUod) 


1     ORIGIN*  TIN  G  ACTIVITY  (Cotporoto  outhor) 

Naval  Postgraduate  School 
Monterey,  California  93940 


2a     REPORT  IICUKITY    CLASSIFICATION 

UNCLASSIFIED  


26  OROUP 


3  REPORT  TITLE 

Application  of  Numerical  Inversion  of  the  Laplace  Transform  to  the  Inverse 
Problem  in  Transient  Heat  Conduction 


4     DESCRIPTIVE  NOTES  (Typo  ot  roport  ond  Incluolvo  dotoo) 

Master's  Thesis 


5    AUTHORfS)  (Loot  nam*,  tint  nam*,  initio!) 

Wendt,  Terrill  J.,  Ensign,  USN 


6     REPO  RT  DATE 

March   1968 


7a.    TOTAL  NO.   OF    PAtll 


119 


7b.    NO.  OF  REFS 


8a      CONTRACT    OR    QRANT    NO. 


b.    PROJECT  NO. 


•  a.    ORISINATOR'S  REPORT  NUMBERfS.) 


tb.   OTHER  REPORT   uo(S)   (A ny  othor  numboro  tttmt  moy  bo  oootfiod 
mio  roport) 


10    AVAILABILITY/LIMITATION  NOTICES 


13     ABSTRACT 

The  direct  problem  in  transient  heat  conduction  requires  the  determination 
of  conditions  at  an  interior  location  when  conditions  are  known  at  the  boundaries 
of  a  solid.   Conversely,  the  inverse  problem  requires  the  determination  of 
conditions  at  the  boundaries  of  a  solid  when  conditions  are  known  at  an  interior 
location.   Consequently  special  methods  are  required  in  the  solution  of  the 
inverse  problem.  A  new  method,  numerical  inversion  of  the  Laplace  transform, 
is  used  to  solve  this  complex  problem.  Application  of  this  numerical  technique 
to  the  semi- infinite  solid,  "long"  cylinder,  and  sphere  is  made,  and  the  accuracy 
of  solution  is  discussed.  This  method  of  solution  provides  the  engineer  with  a 
simple,  powerful  tool  that  can  be  used  in  the  determination  of  heat  transfer 
phenomena  in  a  solid. 


DD  M!*.<  1473 


120 


UNCLASSIFIED 


Security  Classification 


UNCLASSIFIED  - — -__ 

Security  Classification 


KEY  WO  R  OS 


Transient  Heat  Conduction 
Inverse  Problem 
Laplace  Transform 
Numerical  Inversion  Technique 


■ 


-  - 


DD  :™..U73  (back) 


121 


UNCLASSIFIED 


S/N     0101-807-6821 


Security  Classification 


